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ABSTRACT 

Using very deep near-infrared VLT/ISAAC and optical HSTAVFPC2 imaging in the Hubble Deep Field 
South and the field around the cluster MS 1054-03, we study the rest-frame ultraviolet-to-optical colors and 
magnitudes of galaxies to redshift z ^ 3. Whereas the present-day color-magnitude distribution shows a promi- 
nent red and blue sequence of galaxies, we find no evidence for a red sequence at z ^ 3. We do find a well- 
defined color-magnitude relation (CMR) for blue galaxies at all redshifts, such that galaxies more luminous in 
the rest-frame V-band have redder U -V colors. The slope of the blue CMR is independent of redshift with 
an average 5{U -V)/6Mv = -0.09 ±0.01. Spectra and colors of z = comparison galaxies suggest that the 
slope can in principle be explained by a trend of increasing dust-reddening with optical luminosity, with mi- 
nor contributions from age and metallicity. The rest-frame color at fixed luminosity of the blue CMR reddens 
strongly A{U - V) « 0.75 from z ^ 3 to z = 0. Much of the color evolution can be explained by aging of the 
stars, although degeneracies in the models prevent us from disentangling the contributions of age and dust. The 
color scatter of the blue sequence is relatively small criU-V) = 0.25 ±0.03 and constant to z 3. Notably, the 
scatter is asymmetrical with a sharp blue ridge and a wing towards redder colors. We explore a set of ensemble 
star formation histories for blue galaxies to study the constraints placed by the shape of the color scatter at 
z = 2 - 3. Models with purely constant or exponentially declining star formation fail to reproduce the correct 
shape, but episodic star formation - here implemented as a two-state model with high and low star formation 
- reproduces the detailed properties well. The combined constraints on the best-fit parameters of the episodic 
model are that the duty cycle, the fraction of time spent in the active state, is more than 40%, and that the 
contrast in SFR between the high and low phase, i.e. the relative strength of the burst, is more than a factor of 
5, which corresponds to more than 0.35 dex of scatter in \og{SFR) around the mean. Episodic models allow 
blue sequence galaxies to have started forming at z >> 3, relying on recurrent star burst to keep the galaxy 
colors blue and suggesting a well defined blue sequence may be found at higher redshifts. However, episodic 
models do not naturally explain the observed tail of very red galaxies, primarily Distant Red Galaxies (DRGs) 
with observed Js-Ks > 2.3. DRGs may have ceased star formation altogether or are more heavily obscured. 
Finally, the relative number density of red, luminous My < -20.5 galaxies increases by a factor of 6 from 
z = 2.7 to z = 0.5, as does their contribution to the total rest-frame V-band luminosity density. We are likely 
viewing the progressive formation of red, passively evolving galaxies. 

Subject headings: galaxies: evolution — galaxies: high redshift — infrared: galaxies 



1. INTRODUCTION 

In the local universe, the color distribution of galaxies is bi- 
modal, primarily reflecting a relationship between color and 
morphology (e.g., Holmberg 1958, Tully et al. 1982). The red 
peak is populated mainly by spheroidal-like early type galax- 
ies and the blue peak consists generally of disk-like late type 
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galaxies, as has later been verified in detail (Hogg et al. 2003, 
Baldry et al. 2004). 

Both populations obey relations between their integrated 
U -V colors and absolute V-band magnitudes, a red and a 
blue sequence, where the colors become systematically red- 
der with increasing luminosity (Baum 1959; Visvanathan & 
Sandage 1977; Chester & Roberts 1964, Visvanathan 1981; 
Tully, Mould, & Aaronson 1982). The optical colors trace 
the integrated light of the galaxy stellar populations and thus 
any correlation of color with magnitude reflects the trend 
of some combination of star formation history (SFH), initial 
mass function (IMF), metallicity, and dust attenuation with 
optical luminosity. 

The relation for red early type galaxies is well-defined, in 
particular in clusters, and is generally interpreted as a corre- 
lation between metallicity and luminosity (Faber 1973, Lar- 
son 1974). The small intrinsic scatter and slope of this rela- 
tion has been used to place limits on the star formation and 
merging histories of these galaxies (Bower, Lucey, & Ellis 
1992; Schweizer & Seitzer 1992; Bower, Kodama, & Ter- 
levich 1998, van Dokkum et al. 1998) 

The relation for star forming blue galaxies is not as well 
understood as that for passively evolving galaxies on the red 
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sequence. The scatter is larger (Griersmith 1980), its origin 
might be more complex, and the slope has been attributed to 
trends with stellar age (Peletier & de Grijs 1998, Brinchmann 
et al. 2004), dust attenuation (Tully et al. 1998), and/or metal- 
licity (Zai-itsky, Kennicutt, & Huchra 1994, Bell & De Jong 
2000), largely depending on the wavelength studied. 

The constraints on galaxy formation models become in- 
creasingly powerful if color-magnitude trends are studied to 
higher redshift, thus over a long time period. For example, 
Kodama & Arimoto (1997) were able to ruled out correlation 
with age as the main cause for the red CMR (cf. Worthey, 
Trager, & Faber 1995), because the predicted evolution of the 
CM sequence with redshift was more than observed. And Bell 
et al. (2003) defined a photometric red sequence to study its 
evolution since z = 1 finding a x2 build-up of stellar mass, in- 
consistent with a scenario where all early type galaxies form 
at high redshift and passively evolve to the present day. 

In this paper, we will investigate the evolution of the rest- 
frame optical colors and magnitudes to much higher redshift 
z ^ 3. We are interested in quantifying any changes in the 
color-magnitude diagram from z = 3 to z = 1, and in particular 
to see how the red sequence and blue sequence evolve over a 
large fraction of the Hubble-time. 

To study the evolution to z = 1 - 3 at the same rest-frame op- 
tical wavelength requires very deep near-infrared (NIR) imag- 
ing, because the rest-frame V band shifts into the NIR at z > 1 . 
While instruments on large telescopes are now making it pos- 
sible to select fairly large samples of massive galaxies by their 
rest-frame optical light out to z ~ 3 (e.g., McCarthy et al. 
2001, Cimatti et al. 2002, Labbe et al. 2003, Forster Schreiber 
et al. 2006, van Dokkum et al. 2006), very few fields reach 
the depth required to study color-magnitude relations. The 
reason is that the CMR is a relatively subtle feature, and re- 
quires high signal-to-noise photometry over a large range in 
magnitudes. 

We have very deep optical-to-NIR imaging in two 
fields from the F aint Infrared Extragalactic Survey (FIRES; 
iFranx et ani200Q) . an optical-to-infrared multicolor survey of 
NIR-selected galaxies. The two fields are the HDFS, contain- 
ing some of the deepest optical and NIR imaging ever taken, 
and the 5x larger field around the z = 0.83 cluster MS 1054- 
03, reaching « 0.7 mag shallower magnitude limits. The deep 
NIR imaging, the broad wavelength coverage in 7 or 8 filters 
from 0.3/im to 2.2/im, and homogeneous photometry make 
it perfectly suited for studies of rest-frame optical colors and 
magnitudes (see Labbe et al. 2003, Forster Schreiber et al. 
2006). 

We shall focus mainly on the sources that populate the blue 
peak of galaxies to high redshift. Many of these would satisfy 
the so-called "U-dropout" criteria, an effective color selection 
technique pioneered by Steidel et al. (1996a,b). However, 
NIR-selection provides a more complete census at z = 2 - 3 
that also includes the reddest galaxies, such as the Distant 
Red Galaxies (Franx et al. 2003, see also Daddi et al. 2004), 
ensuring that any color trend we find is not a selection ef- 
fect. The ultradeep surveys in the Hubble Deep Field North 
(HDFN) and Hubble Deep Field South (HDFS) have already 
shown that z = 2-3 blue galaxies obey a clear color-magnitude 
relation, qualitativel y similar to the relation for blue late-type 
in the local univer se dPapovich. Dickinson. & Fergusonl2001l : 
iLabbe et"ani2003h . We will attempt to establish a direct con- 
nection, and investigate which other properties of the relation, 
such as the scatter, can be used to constrain the star formation 
histories of blue galaxies to z = 3 and higher. 



This paper is organized as follows. We present the data 
in §2, describe in §3 the general rest-frame color-magnitude 
distribution of galaxies up to redshift z ^ 3, we quantify the 
evolution of the blue color-magnitude relation in §4 and study 
its origin in §5. We use the observed scatter of blue sequence 
galaxies to place constraints on models for galaxy evolution 
in §6, and quantify the onset of the red galaxies since z = 3 
in §7. Where necessary, we adopt an Sl^ = 0.3, J^a = 0.7, 
and Ho = 70 km s"'Mpc"' cosmology. We use magnitudes 
calibrated to models for Vega throughout. 

2. THE DATA 
2.1. The Observations and Sample Selection 

The observations were obtained as p art of the public F aint 
Infrared Extragalactic Survey (FIRES; Franx et al . 2000) the 
deepest groundbased NIR survey to date. It covers two fields 
with existing deep optical WFPC2 imaging from the Hub- 
ble Space Telescope (HST): the WPFC2-field of HDFS, and 
the field around the z = 0.83 cluster MS 1054-03. The obser- 
vations, data reduction, and assembly of the catalog source 
catalogs ar e presented in detail bv Labbe et al.l (l2003h for the 
HDFS and lForster Schreiber et alJ (l2006h for the MS 1054-03 
field. We will summarize the main steps here. 

The two fields were observed in the NIR /,,//, and bands 
with th e Infrared Spectrometer and Array Camera (ISAAC; 
iMoorwood 1997) at the Very Large Telescope (VLT). In the 
HDFS we spent a total of 101.5 hours in a single 2.5' x 2.5' 
pointing, resulti ng in the dee pest public groundbased NIR 
imaging to date jLabbe et al.l 12003). We complemented the 
data set with existing ultradeep optical imaging from HST in 
the 1/300,8450, ¥(,()(,, hi4 bands (Casertano et al. 2000), where 
the passbands refer to the HST/WPFC2 F300W, F450W, 
F606W, and F814W filters, respectively. In the field of 
MS 1054-03 a further 77 hours of NIR imaging was spent on 
a mosaic of four ISAAC pointings centered on the z = 0.83 
cluster MS 1054-03 (Forster S chreiber et al. 2006). We com- 
plemented these data with WFPC2 mosaics in the V(,o6 and 
/814 bands ( van Dokkum et all2000h . and collected additional 
imaging with the VLT FORSl instrument in the U,B, and V 
bands. In both surveyed fields the seeing in the final NIR im- 
ages was 0"45-0f'55 FWHM. 

In each field, we registered the images and convolved them 
to match the image quality of the band with the worst see- 
ing. We detect objects in the Ks-hand using versi on 2.2.2 of 
the SExtractor software ((Ber tin & Arnoutslri996l) . We per- 
formed photometry in matched apertures, measuring colors 
in customized isophotal apertures defined from the /Q-band 
image, and measuring total ^T^-band magnitudes in a cus- 
tom ized elliptical Kro n ap erture (SExtractor MAG AU TO). 
See'Labbe etal] (l2003h and lForster Schreiber et all ( 12006 ) for 
details. Random photometric uncertainties were estimated 
empirically from the flux scatter in apertures placed randomly 
on empty parts of the map. Systematic calibration uncertain- 
ties between filters are estimated to be less than 5%. The 
total 5-a limiting depth for point sources is K'"'" = 23.8 for 
the HDFS, and 23.1 for the MS1054-field. The source cata- 
log contains a total of 833 sources from the HDFS and 1858 
sources from the field of MS 1054-03. 

2.2. Redshifts and Rest-Frame colors 

For most sources we rely on photometric redshifts as spec- 
troscopic redshifts are available for a small subset only. Pho- 
tometric redshifts were estimated by fitting a linear combina- 
tion of redshifted empirical and model galaxy spectra to the 
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observed flux points. The template set ranges from very blue 
10 Myr (iBruzual & Charlo j l2003l) burst models to very red 
empirical elliptical templ ates. The a lgor ithm is described i n 
detail by Rudnick et al. (2 00Tll2003h and lLabbe etall (l2003h . 
We only determine photometric redshifts for sources with the 
highest-quality photometric data. Our quality cut requires that 
a source must have photometric information in all bands and 
a minimal nominal exposure time of 20%. We adopt a min- 
imux flux error of 5% for all bands to account for absolute 
calibration uncertainties and for mismatches between the ob- 
servations and the adopted galaxy template set. 

We derived photometric redshifts for the 1475 out of 2691 
sources that met our quality criteria. Since redshift errors are 
a primary source of uncertainty in the rest-frame luminosities 
and colors they must be well charactarized. To estimate pho- 
tometric redshift errors we performed Monte-Carlo simula- 
tions, randomizing the observed fluxes within their errors and 
finding the best-fit redshift again. The method accounts for 
the effect of photometric uncertainties, template mi smatch, 
and the possibility of secondary solutions (see IRudnick et alJ 
|2003). For sources with both good photometry and spectro- 
scopic redshifts we find good agreement between the Monte- 
Carlo redshift error Szph.MC and the difference between photo- 
metric and spectroscopic redshift Szspec =< \zspec - Zphot 1/(1 + 
Zspec) >■ The mean photometric redshift uncertainty for all 
/Tj-selected from direct comparison to available spectroscopy, 
is Szspec = 0.07, while the accuracy for sources at Zspec > 2 is 
better Szspec = 0.05. We identified and rem oved stars using the 
method described in lRudnick et al.l (l2003h . 

We derive rest-frame luminosties U^"' from the observed 
SEDs and redshift information. We estimate U^"' by inter- 
polating between the observed fluxes, using the best-fit tem- 
plates as a guid e. Details on deriving L"" are described 
extensively by (iRudnick et al.l l2003h . As our rest-frame 
photometric filter system we use the ultraviolet HST/FOC 
F140W,F170W, and F2 20W filters and the optical filters 
UX,B and V of |Bessel (Il990l) . which sh aU be denoted as 
1400,1700,2200,f/,B, andV. The|Besse| system was 

calibrated to the Dreiling and Bell (1980) model spectrum for 
Vega, while the HST/FOC system was calibrated to the Ku- 
rucz (1992) model for Vega. 

The rest-frame luminosities and colors are sensitive to the 
uncertainties in the photometric redshifts. Therefore, in the 
remainder of the paper we only analyze the sample of galax- 
ies with spectroscopic redshifts or those with reasonable red- 
shift uncertainties 5zph.MC /i^+Zph) < 0.2, keeping 1354 out of 
1475 galaxies. Potentially, such a quality cut could introduce 
a bias, e.g., if blue galaxies were to have systematically dif- 
ferent photo-z undertainties than red galaxies. However, we 
find no such color-dependence for our sample and we verified 
that the rest-frame color distribution of the rejected galaxies 
is consistent with that of the galaxies we kept in the sam- 
ple. Because our analysis is sensitive only to the rest-frame 
color distribution we expect no significant bias. The median 
Szph,Mc/i^+Zph) for the remaining galaxies is 0.05. The re- 
duced images, photometric catalog, redshifts, and rest-frame 
luminosities are all available on-line through the FIRES web- 
site". 



3. THE REST-FRAME COLOR-MAGNITUDE DISTRIBUTION OF 
GALAXIES TO Z ~ 3 

' ' http://www.strw.leidenuniv.nl/~fires 



The ultraviolet-to-optical color-magnitude distribution is a 
well-studied diagnostic in studies of low redshift galaxies 
(e.g., Strateva et al. 2001; Hogg et al. 2002; Baldry et al. 
2004). Quantifying this distribution out to higher redshifts 
will provide strong constraints on models of galaxy forma- 
tion, which must reproduce these observations (e.g.. Bell et 
al. 2004, Giallongo et al. 2005). 

In Fig. 1 we present the rest-frame 2200 -V colors versus 
absolute V magnitude of galaxies in the FIRES fields. Here 
we use the 2200-band instead of the more common [/-band 
as it is traced by our deepest HST observations to z ^ 3. To 
show the evolution with time we divided the galaxies in three 
redshift bins centered on z = 2.7, 1.8, and 1.1, which corre- 
sponds roughly to 2.4, 3.6, and 5.6 Gyr after the big bang with 
a time-span of 1, 1.5, and 2.7 Gyrs, respectively. The redshift 
ranges are defined so that the rest-frame 2200 and V-band lie 
in between our observed filters, while the widths of the bins 
are a trade-off between reducing the effects of evolution over 
the time-span of the bins, while keeping a statistically signifi- 
cant sample in each interval. 

Obvious from Fig. 1 is that at all redshifts galaxies occupy a 
fairly narrow blue locus in color-magnitude space. The z < 1 
blue locus consists of late-type spiral- and irregular galaxies, 
whi le the z > 2 blue locus is populated by tZ-drop out galax- 
ies (iMadau et al.lll996l;lGiavalisco & Dickinsonll20bli Steidel 
et al.l996a,b)''^. The blue peak has a well-defined ridge on 
the blue side, while galaxy colors spread out several magni- 
tudes to the red. Furthermore, the blue ridge is tilted: more 
luminous galaxies along the blue ridge tend to have redder 
2200 -V colors. This ridge is what we define as the color 
magnitude relation of blue galaxies. At the brightest magni- 
tudes there is some evidence for an upturn of the blue CMR, 
where most galaxies lie to the red of the relation defined by 
faint galaxies. Baldry et al. (2004) found qualitatively similar 
results at z ~ 0.1 for galaxies in the Sloan survey. The knee 
of this upturn evolves from roughly My ^ -22.5 at redshift 
z 3 to -2 1 at z 1 . The scatter around the blue CMR 
is markedly asymmetric, as can be clearly seen at z ~ 3. The 
scatter has a blue ridge and a red wing, with which we mean 
the skew towards red colors near the peak of the distribution. 

It is unlikely that the blue CMR is artificial. One obvious 
worry is that if the templates used in the photometric redshift 
code are all redder than the true colors of the galaxies, then the 
photometric redshift algorithm could create an artificial color 
magnitude relation with a blue ridge. However, our template 
set reaches much bluer colors than what we find for the galax- 
ies. Magnitude selection effects are also unlikely to cause or 
affect the blue CMR; we are not biased against galaxies that 
are bright in My and blue in 2200 - V and our photometry 
is deep enough to ensure the trend is not caused by an ap- 
parent lack of faint red galaxies. At low redshifts z < 1 we 
may perhaps start to miss some the faintest blue galaxies, as 
the Ks-hand selection probes significantly redder wavelengths 
than the redshifted V-band, but this effect is small and does 
not affect the analysis presented in this paper. 

The color magnitude distribution has an extended tail to 
very red colors at all redshifts. The red tail extends up to 4 

(/—dropout gal axies in the HDFS field were selected with the HST cri- 
teria offaiavalisco & Dickins og ( 2001) and (/-dropouts in the MS 1054-03 
field were selected with(;7-B)> IAB-V)< 1.4, and ((/-B) >(B-V)+0.5. 
This is not a proper U-dropout selection per se, as the VLT/FORS UBV fil- 
ter set does not map directly to the Palomar UnGR filterset, but galaxies in 
our sample with these UBV colors have synthe tic U nGR colors wh ich would 
satisfy the conservative "C" and "D" criteria of lSteidel et alj j20O3l) 
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Fig. 1 . — The rest-frame 2200- V colors versus absolute V— band magnitude for galaxies in the field of the HDFS (left) and in the field of MS 1054-03 (right). 
The observations are split into three redshift bins with the highest redshift on top. The errorbars denote the Itr uncertainties on the rest-frame colors. Galaxies 
redward of the dotted line have uncertain rest-frame colors, but can still be observed. Distant Red Galaxies (Franx et al. 2003), selected by their red observed 
Js-Ks > 2.3 colors are mai'ked in red, [/-dropout galaxies are shown in blue. The solid line shows a fit of a linear relation with a fixed slope of -0.17 to the 
galaxies in the blue peak of the color-magnitude distribution. The thick gray fine marks the My Hmit corresponding to our ^C,— band magnitude limit, calculated 
by adding to the My of each source the quantity (Ksource — Ku,,,!, ) hence taking into account the color dependence of the rest-frame detection Hmit. The inset 
panels show histograms of the residual color distribution after the slope is subtracted (white histograms). The peak of the distribution evolves to redder colors 
from high to low redshift. The bright end of the galaxy distribution, to My < -19.5 (HDFS) and My < —20.5 (MS1054), evolves rapidly in color between z ~ 3 
and z ~ 1 (gray histograms). The tentative build-up of a second, red peak of galaxies can be seen as early as z = 1.4— 2. We note that the field of MS 1054 contains 
a massive cluster atz = 0.83 which significantly enhances the number of bright red galaxies in the z ~ 1 bin. 



magnitudes and is bound in the red at low redshift by a red se- 
quence of galaxies. The onset of the red sequence can be ten- 
tatively observed as early as z 2 in the field of MS 1054-03 
and becomes clearly visible in in the MS 1054-03 field at z ~ 
1 . The narrow red sequence here is largel y attributable to the 
elliptical galaxies in the cluster at z = 0.83 (lvan Dokkum et al.l 
HOOO). 

Finally, the color histograms in Fig. 1 show that at the bright 
end the galaxy population undergoes a strong evolution in 
color as the relative number of luminous blue galaxies de- 
creases with time and that of red galaxies increases. As our 
galaxy sample is obviously too small to describe properly the 
upturn at the bright magnitudes of the blue sequence or the 



evolution of the red sequence, the primary focus of our anal- 
ysis will be on the linear part of the blue CMR. 

4. THE COLOR-MAGNITUDE RELATION OF BLUE FIELD 
GALAXIES 

We now quantify the three main properties of the blue 
CMR: the slope, the zeropoint, and the color scatter around 
the relation. In addition, we describe the evolution with red- 
shift. 

4.1. The Evolution of the Blue Sequence Slope 

While the distribution of UV-to-optical colors and absolute 
V-band magnitudes of faint blue galaxies may be well char- 
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TABLE 1 

The Blue Color Magnitude Relation 







HDFS 




MS 1054 








slope'' 


intercept 


slope'' 


intercept'^ 


0.5- 


0.7 


-0.176 


0.226 


-0.262 


0.098 


0.7- 


1.4 


-0.206 


-0.186 


-0.144 


-0.190 


1.4- 


2.2 


-0.140 


-0.705 


-0.234 


-0.516 


2.2- 


3.2 


-0.139 


-0.950 


-0.324 


-0.892 



Note. — 
" The redshift range of the subsample 
* The best-fit slope of the 2200 -V versus V-band CMR. 
^ The intercept of the blue CMR at My = -21, fit with a slope fixed at - 



-0.17. 



Fig. 2. — The evolution of the slopes of the blue color-magnitude rela- 
tion (CMR), derived from linear fits to the color-magnitude distribution in 
the HDFS {filled circles) and MS 1054-03 field (diamonds). The slopes are 
plotted at the mean redshift of the galaxies in the bin. The uncertainties de- 
note the 68% confidence interval obtained from bootstrap resampling. The 
2200 — y slope is constant up to z ~ 3. We show the predicted slopes for 
blue sequence galaxies in the Nearby Field Galaxy Survey (NFGS; Jansen, 
Franx, & Fabiicant 2000a) if the local blue CMR were caused by dust red- 
dening (star) or stellar age (triangle). Also shown is the predicted slope of 
the 2200 — V slope for red, early-type galaxies in the Coma cluster assuming 
a metallicity-luminosity relation (Bower, Lucey, Sl Ellis 1992). 



acterized by a linear relation, this is not true for all galax- 
ies, certainly not at lower redshifts. As a result, a straight- 
forward least-squares linear fit to all galaxies will fail to esti- 
mate the blue slope, unless red outliers are removed from the 
fit. "Robust" techniques, such as iteratively rejecting outliers 
(sigma clipping), or assigning low weights to outliers (e.g., 
the biweight estimator of Beers, Flynn, and Gebhardt 1990) 
do not perform very well either because the relative number 
red galaxies changes significantly with redshift and becomes 
fairly large at z < 1 . 

Therefore, we used an even more robust technique, called 
"mode-regression", which works as follows. For a range of 
values for the slope, we calculate the distribution of residual 
colors relative to that slope. We choose as our best fit slope the 
value that results in the highest, narrowest peak in the residual 
color-distribution. The peak height is calculated using a gaus- 
sian kernel density estimator; a smoothed histogram where 
each data point is replaced by a gaussian kernel. The width 
of the kernel is set to the median uncertainty in the rest-frame 
galaxy colors. The mode-regression technique is very insen- 
sitive to outliers as long as most galaxies are on the blue se- 
quence, a condition that our sample satisfies. There is a small 
systematic dependence on the adopted width of the smooth- 
ing parameter. A large smoothing paramater flattens the slope, 
but in our case this systematic effect is small compared to the 
random uncertainties. We estimated uncertainties in the fit 
parameters by bootstrap resampling the color-magnitude dis- 
tributions 200 times with replacement, repeating the fitting 
procedure, and taking the central 68% of the best-fit parame- 
ters as our confidence interval. 

Figure|2]shows the evolution of the slopes 5{22Q0-V) /SMy 
with redshift (see also Table 1). The 2200 -V color is used, 
rather than U -V , because it has higher signal-to-noise and a 
larger wavelength baseline, which helps to detect evolution of 
the slope. The U -V slope is also given to enable direct com- 
parison to other work. We plot the results in the fields of the 
HDFS and MS 1054-03 separately, to demonstrate the results 



are consistent in the independent fields. The measured evolu- 
tion of the slopes from all our measurements is 0.00(±0.01)z 
in U-V or 0.01(±0.03)z in 2200-y over 0.5 < z < 3, con- 
sistent with zero. The error-weigthed mean is: 

(5(t/-y)/5My = -0.09±0.014 (1) 
(5(2200-y)/5Mv =-0.17 ±0.021 (2) 

We compare the U-V slope to local observations of blue 
galaxies from the Nearby Field Galaxy Survey (NGFS; Jansen 
et al. 2000a). Using the same technique as described above 
we find (5(f/-y)/(5Mv =-0.078 ±0.014, consistent with the 
high-redshift measurements. Coincidently, the U-V slope 
of nearby early-type galaxies in the Coma cluster 5{U - 
y)/(5Mv =-0.08 ±0.01 also has a similar value. 

To compare the local results to our 2200-y measurements 
we transform the local U magnitude into 2200 magnitudes. 
As expected, the exact transformation depends on the assump- 
tion about the origin of the slope. If the slope is caused 
by dust flie tr ansformation is A(2200-y) = 2.28A(f/-y) 
(ICalz etti et al.l 120001 reddening). In case of age, with more 
luminous galaxies having older stellar po pulations, the trans- 
formation is A(2200-y) = 1.6A(f/-y) (iBruzual & Chariot! 
i2003i models). This transformation holds to better than 8% 
for blue stellar continuum colors {U -V) < 0.5 over a range 
of declining star formation histories (e-folding timescales r > 
1 Gyr) and metalUcities (Z = 0.02-0.008). Note, in the pres- 
ence of any dust reddening, the intrinsic stellar continua are 
even bluer, and the approximation would be even more accu- 
rate. The thus transformations result in S(220Q-V)d„it /SMy = 
-0.17±0.05, and 5(2200-y)flg,/(5Mi/ = -0.12±0.03. The 
red galaxy color-magnitude relation from the Coma cluster 
can be understood as a correlation between luminosity and 
metallicity. Fitting old passively evolving BC03 model spec- 
tra and a range of metalUcities, we find 6(2200 -V)/SMv = 
-0.08 ±0.01. 

We conclude that we find no evidence in our data for evo- 
lution of the U-V and 2200 - V slope from redshift z = to 
z ^ 3. Only if we use the 2200-y slope inferred for red early- 
type galaxies in Coma is there a hint of evolution at z < 1 . This 
is implausible however, as early type galaxies have much red- 
der colors and higher stellar ages than blue sequence galaxies, 
and this option is not considered further. 

4.2. The Blue Sequence Slope at z=3 as a Function of 
Rest-Frame Color 

There is no need to limit the analysis to only U -V or 
2200 - y colors as the dataset probes a much wider range 
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Fig. 3. — The blue CMR slope at z = 2 — 3 versus rest-frame wavelength 
in the field of the HDFS. The points show the observed slope in the rest- 
frame A - y color versus My as a function of the filter A, where A is the 
1400, 1700, 2200, (/, and B-band. Overplotted are expectations for three red- 
dening laws scaled to fit the U,B, and V points: the Calzetti et al. (2000; 
COO) dust law (dashed line), the MW extinction law (Allen 1976; dotted), 
and the SMC extinction law (Gordon et al. (2003); dash-dot line). The thick 
solid Une shows the color dependence in the case that stellar population age 
correlates with My (Solar metallicity, constant star formation). The gray area 
shows exponentially declining SFHs with r > 1 Gyr. 



Fig. 4. — The evolution with redshift of the zeropoint of the blue CMR at 
a fixed My = — 21. Symbols are as in Fig. 1. The lines represent tracks of 
Bruzual & Chariot (2003) stellar population models at a fixed My; hence the 
models have been corrected for luminosity evolution. We show a model with 
formation redshift Zf = 3.2, a star formation timescale t = 10 Gyr, and fixed 
ICalzetti et all gOOO) reddening of E(B-V) = 0.15 (dashed line); one with 
Zf = 10, constant star formation, and fixed E(B — V) = 0.15 (solid line). We 
also show a model with zf = 10, r = 30 Gyr, and E(B — V) evolving linearly 
in time from at z = 10 to 0. 15 at z = (dotted line). While most of the color 
evolution can be attributed to aging of the stellar population, the simplest 
star histories do not fit the color evolution very well. Better agreement can 
obviously be obtained with more complex SFHs or with variable levels of 
dust reddening. 
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in rest-frame wavelengths. Hence we investigate the depen- 
dence of the blue CMR slope on rest-frame color. Of partic- 
ular interest are the wavelengths spanning the age-sensitive 
3650 A (Balmer) and 4000 A breaks, and the wavelengths 
covering the far-UV, where the slope is sensitive to dust red- 
dening. Therefore, the following discussion is limited to the 
HDFS dataset, which is the deepest, and to only the highest 
redshift bin (2.2 < z < 3.2), where our observed filter set di- 
rectly probes the wavelengths of interest (1400-5500 A). 

We calculated the rest-frame A - V colors, where A is 
the 1400, 1700, 2200, t/, or B-band (see §2.2), and fitted the 
CMR slope for each color versus My as in §4.1. Figure 3 
shows the blue slope as a function of color. Clearly, the 
A- V slope changes as A mo ves to bluer wavelengths. Us- 
ing lBruzual & Chariot! ( l2003h models, we overplot in gray the 
predicted color dependence if the U-V slope were caused by 
trends with age or star formation history. In addition, we show 
several common reddening curves scaled with arbitrary con- 
stants to fit the f/-y slope {5E{B-V)/5Mv = -0.04 for a 
Calzetti et al. 2000 law; -0.05 for the Milky Way extinction 
curve, Allen 1976; and -0.02 for the SMC curve, Gordon et 
al. 2003). 

As can be seen, an age or SFH related slope has a differ- 
ent wavelength dependence than the data. This model does 
not fit at all to the bluest point, which is the CMR slope in 
1400-y color versus My- In contrast, the predicted Calzetti 
et al. (2000) and SMC reddening curves lie very close the 
data, which may indicate that the slope is primarily caused 
by a trend with dust. It may also imply that the Calzetti 
reddening law is appropriate for z ~ 3 galaxies, at least for 
galaxies on the blue sequence. The MW curve fits worse as 
the characteristic bump at 2175 A is not seen in the data. 
Mixed dust geometries can probably dilute the bump, al- 
though Gordon, Calzetti, & Witt (1997) suggest the lack of 
the bump feature in nearby starburst galaxies is probably in- 



trinsic, perhaps due to the high UV energy densities in the 
star forming regions. In the case of Calzetti or SMC red- 
dening, the implied systematic change of Ay along the blue 
CMR from My = -19 to My = -23 is < Ay -Ay,,„i„ < 0.7 
and < Ay-Ay,,yii„ < 0.3, respectively, where the extinction at 
My = -19 is Ayjnin- We notc that the distribution in Ay values 
inferred from SED fitting of individual galaxies can be sub- 
stantially br oader due the degeneracy between age and dust in 
the models (Pap ovich. Dickinson. & Ferguson ^2001, e.g.,), or 
sample selection, e.g., the sample of tZ-dropout galaxies of 
Shapley et al. (2001) is biased to galaxies with larger implied 
dust obscuration. 

4.3. The Evolution of the Blue Sequence Zeropoint 

As we find no evidence that the slope of the blue se- 
quence evolves with redshift, we fix from hereon the slopes 
at 5{U-V)/5Mv = -0.09 and 6{220Q-V) / 5My = -0.17 as 
determined from our data. To find the zeropoint of the rela- 
tion, we subtract the linear relation shifted to My = -21 and 
define blue sequence zeropoint as the location of the peak of 
the residual color distribution. The 68% confidence interval 
of the zeropoint is obtained with bootstrapping. We interpret 
the zeropoint as the average color of galaxies on the blue se- 
quence with absolute magnitude Mv/ = -21. 

Figure |4] presents the evolution of the 2200 - V zeropoint 
with redshift. We also plot the inferred 2200- V color of local 
blue sequence galaxies from the NFGS, using BC03 models 
wi th a constant star formation history and E(B-V) = 0.13 
of ICalzetti et alj (1200 0) reddening. The reddening is ap- 
propriate for nearby My = -21 blue seque nce galaxies (see 
iJansen. Franx. Fabricant. & Caldwellll2000l and §5.1). Again 
we will also give the evolution of the U-V zeropoint. Clearly, 
the color of the blue CMR at fixed absolute magnitude red- 
dens monotonically from z ^ 3 to z ~ 0.5. The galaxies be- 
come redder in fZ-V by « 0.5 mag from z ~ 2.7 to z ^ 0.5, 
and in 2200 -V by w 1.1 mag (from the FIRES data alone). 
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TABLE 2 
The Scatter around the Blue CM 







HDFS 




MS 1054 












^obbs 






0.5-0.7 


0.47 


0.45 


0.56 


0.54 




0.7-1.4 


0.58 


0.56 


0.59 


0.57 




1.4-2.2 


0.40 


0.37 


0.48 


0.46 




2.2-3.2 


0.51 


0.49 


0.62 


0.60 



Note. — 
" The redshift range of the subsample. 

* The biweight scatter in 2200 — V color after removing galaxies that are 
A(2200-y) > 1.6 redder than the CMR. 
The intrinsic color distribution after subtracting the observational errors. 



Fig . 5 . — The observed scatter of galaxies on the blue CMR as a function of 
redshift. Symbols as in Fig. 1 . There is no evidence for evolution of the scat- 
ter with redshift. The dashed Une denotes the mean of all measurements. The 
gray dotted line denotes the median photometric error of individual galaxies. 
The photometric errors are small in comparison, showing that the observed 
scatter is intrinsic, not caused by photometric errors. 



The total variation including the z = point is 5{U-V) ~ 0.75 
and (5(2200 -y)« 1.4. 

Surprisingly, a straight hne describes the points at z = 0.5-3 
rather well and we find 

f/-y = 0.45(±0.06)-0.29(±0.04)z (3) 
2200 - y = 0.34(±0. 1 ) - 0.50(±0.06)z (4) 

The colors in the independent fields of the HDFS and 
MS 1054 agree very well, suggesting that field to field vari- 
ations do not affect this result. If the empirical linear relation 
derived from our high-redshift data is extrapolated to z = 0, it 
coincides with the z = point measurement from the NFGS. 
This is encouraging because absolute calibration between dif- 
ferent surveys is rather difficult, and the NFGS data has been 
transformed to 2200 -y colors from other passbands. 

We can compare the observed color evolution directly to 
simplistic predictions from stellar populations models if we 
assume that the galaxies remain on the ridge of the CMR 
throughout their life. This assumption may very well be 
wrong, but it enables us to get a feel what range in colors 
basic stellar populations span. In §5.1 and §6 we will attempt 
more elaborate modeling, where we add newly formed galax- 
ies to the blue CMR, and let older galaxies evolve away from 
the relation at later times. Obviously for such simple models, 
the galaxies all have the same color, and these follow directly 
from the colors of the stellar population model, depending 
on star formation history and the dust absorption onl y. Us- 
ing colors and luminos ities from [Bruzual & Chariot ( l2003h 
models and a fixed ICalzetti et al.l (^000) dust reddening, we 
generate colors at a constant absolute magnitude by apply- 
ing a small color correction. It reflects the fact that galaxies 
that were brighter in the past populated a different part of the 
CMR. Using the measured blue CMR slope to apply the cor- 
rection, we find the total amplitude of the effect is less than 
A(t/-y) < 0.05 mag and A(2200-y) < 0.1 mag for most 
models. 

The tracks in Figure|4]indicate that the color evolution from 
z ^ 3 to z can in principle be fully attributed to aging of the 
stellar populations without any dust evolution. However, the 
detailed shape of the evolution is not well reproduced. Models 
with constant star formation fit remarkably badly; the evolu- 



tion is much slower than observed. Models with exponen- 
tially declining SFRs (timescale r = 10 Gyr, formation red- 
shift z/ = 3.2) agree better, but we find no satisfying fit. The 
slope is either too flat (for r > 10) or too red (for r < 10). 
Note that we must restrict the fit to zf > 3.2 as the color- 
magnitude relation is already in place at zf < 3.2. 

Naturally, the models can be made to fit better by allowing 
more freedom to the parameterization of the star formation 
rate or allowing varying amounts of dust. For example, the 
dotted line in Figure|4]shows a r = 30 Gyr model with forma- 
tion redshift z/ =10 and with a reddening evolving linearly in 
timefrom£'(B-y) = atz/= 10to£'(B-y) = 0.13 atz = 0. 
The conclusion is that the amount of color evolution in the 
interval < z < 3 can in principle be ascribed to aging of the 
stellar populations, but there is no unique answer for the aver- 
age star formation hi story, due to the degeneracy age and du st 
in the models (e.g.. lPapovich. Dickinson. & Fergusonll2001h . 

4.4. The Evolution of the Blue Sequence Scatter 

As can be seen in Figure 1, both FIRES fields show a large 
range in galaxy colors up to z 3, more than 4 mag in the rest- 
frame 2200- y color. Yet, the majority of galaxies populate a 
blue color-magnitude sequence with remarkable low scatter 

To measure the scatter around the CMR we again subtract 
from the observations the linear fits obtained with the fixed 
slope (see §4.2 and §4.3). We removed galaxies that are red- 
der than A(t/-y) > 0.5 mag or A(2200-y) > 1.6 mag rel- 
ative to the blue CMR to prevent the build-up of red galaxies 
with redshift to bias the scatter esti mates. As a scatter estima- 
tor w e use the biweight estimator jBeers. Flvnn. & Gebhardtl 
IT990h . which assigns lower weight to points that are far from 
the center of the distribution. For Gaussian distributions this 
estimator reduces to the conventional standard deviation. 

Figure |5] shows the 2200-y scatter (see also Table 2). The 
scatter of the colors around the blue CMR is relatively nar- 
row and constant with redshift up to z ^ 3 with values of 
f7i,v.(2200-y) = 0.51 ±0.05 mag and atAU-V) = 0.25 ±0.03 
mag. Note that the range of the color distribution of all galax- 
ies is much larger (~ 4 mag in 2200-y, see Fig. 1). 

The scatter includes a small contribution of observational 
errors. The median photometric error of the galaxies is 0.10 
and0.14mag in U-V and 2200-y respectively, shown by the 
gray dotted line in Fig.|5] Using only galaxies with observa- 
tional errors smaller than 0.10 and 0.14 in U-V and 2200-y, 
respectively, we find the same scatter. We conclude that the 
measured scatter is predominantly intrinsic, likely the result 
of variations in dust and stellar popultion ages of the galaxies. 
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Fig. 6. — The color scatter perpendicular' to the blue CMR at z = 2-3 
versus rest-frame wavelength in the field of the HDFS. The points show the 
observed scatter in the rest-frame \ -V color as a function of the filter A, 
where A is the 1400, 1700,2200,1/, and B-band. Overplotted are expecta- 
tions for three reddening laws scaled to fit the U,B, and V points assuming 
the scatter is produced by dust alone: the Calzetti et al. (2000; COO) dust law 
(dashed line), the MW extinction law (Allen 1976; dotted), and the SMC ex- 
tinction law (Gordon et al. (2003); dash-dot line). The thick solid line shows 
the color-dependence in the case that stellar population age causes the scatter. 
The gray area shows exponentially declining SFHs with t > 1 Gyr. 



4.5. The Blue Sequence Scatter at z=3 as a Function of 
Rest-Frame Color 

There is some evidence from the broadband colors that the 
scatter at z ^ 3 might be mostly due to age variations. To illus- 
trate this point, figure |6] shows the wavelength dependence of 
the scatter. We first plot the predictions for common dust laws, 
where the curves, scaled to fit the U,B,V data, show the pre- 
dicted FUV scatter if the U-V scatter were purely a variation 
in attenuation. However, the observed (A-V) scatter increases 
much less with decreasing wavelength; the scatter in FUV 
colors is rather small. The wavelength dependence is much 
closer to predictions of age variations, here calculated from 
BC03 stellar populations models for a range of star formation 
histories with declining star formation rates (r > IGyr"'). 

Another way of showing the different behaviour of the 
blue CMR slope and scatter with wavelength is to plot the 
1400-2200 and U-V color magnitude relations side-to-side 
(Fig-IZJ- The scatter (relative to the slope) is 3 times smaller in 
1400-2200 than in the U-V. If the blue CMR slope and scat- 
ter were both caused by the same process, for example dust 
variations, then this ratio should be wavelength independent. 
An explanation is that the U-V scatter is predominantly an 
age effect: U-V colors of stellar populations models evolve 
strongly with age due to the onset of the Balmer/4000A break, 
but less so in FUV colors. 

We can use the scatter in the FUV to place a simple upper 
limit to the contribution of dust extinction to the U-V scatter. 
The observed FUV scatter (!;„,( 1400 -2200) = 0.19 is only 
slightly larger than the estimated photometric uncertainties. 
Assuming it is entirely i ntrinsic this translates to (t(U -V) = 
0.15 for a Calzetti et aP (2000) dust-law, or almost half the 
observed value. 

We will discuss implications of the constant scatter with 
redshift in §5.3, and discuss detailed modeling of its shape at 



•i^ 

o 
o 

CJ 
CJ 

I 

o 

^ -2 




> 
I 

D 



(a) 



'I I ' 



■(b) 



|| I J h 




-24 -23 -22 -21 -20 
My (rest-frame) 



19 



Fig. 7.— The FUV 1400-2200 versus Mv relation (a) compared to the 
U — V versus My relation (fo) at z = 2 — 3 in the HDFS field. Open points 
indicate sources with uncertain 1400-2200 colors. The dashed line in (a) is 
the predicted 1400 — 2200 slope if an age-My correlation (with r = 1 Gyr 
SFH) is fitted to the U — V points, underpredicting the observed slope (see 
Fig 3). Note that the scatter (relative to the slope) is 3 times smaller in 1400- 
2200 color than inU — V, whereas this ratio is expected to be independent of 
wavelength if both slope and scatter are caused be the same process (compare 
Fig 3 and Fig. 6). 



z ~ 3 in §6. 

5. THE ORIGIN OF THE BLUE SEQUENCE 

The cause of the blue color-magnitude relation is impossi- 
ble to determine from the broadband photometry of individ- 
ual galaxies alone. Fits of models to the photometry produce 
stellar age and dust estimates which are very uncertain (e.g., 
Shapley et al. 2001, Papovich et al. 2001). Unfortunately, 
more direct spectroscopic age and extinction determinations 
are notoriously difficult for distant galaxies. The rest-frame 
optical spectroscopy is generally not deep enough for the de- 
tection of H/3 or for accurate measurements of Ha equivalents 
widths over a large magnitudes range. In addition, the blue 
CMR is a relatively subtle feature, hence we can only inves- 
tigate it properly for nearby galaxies that have high signal-to- 
noise integrated spectra and photometry over a large magni- 
tude range. 

5.1. The Origin of the Blue Sequence in the Local Universe 

We analyzed the spectra of a sample of 196 nearby galax- 
ies from the NGFS dJansen et al.ll2000l) . a spectrophotometric 
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Fig. 8. — U-V colors versus absolute V-band magnitude of nearby normal 
galaxies from the Nearby Field Galaxy Survey. Insets are as in Fig. 1. The 
solid line denotes a linear fit to the blue sequence {solid line). We select 
blue sequence galaxies by requiring colors within A{U-V) < 0.5 mag of the 
linear fit (dashed line). 



survey of a subsample of galaxies from the CfA redshift sur- 
vey (Huchra et al. 1983). The sample matches the distribu- 
tions of morphology and magnitude of the nearby galaxy pop- 
ulation, and covers a large magnitude range -15 < Mb < -23. 
The NFGS provides integrated broadband photometry and 
integrated spectrophotometry of all galaxies, including line 
fluxes and equivalent widths of [Nil], Hq, and Hp. 

We chose not to use the Sloan Digital Sky Survey as our 
low redshift comparison as the SDSS spectroscopic sample 
covers only a small observed magnitude range (« 2.5 mag) 
in a fixed small aperture (Strauss et al. 2002). Delicate com- 
pleteness and aperture c orrections (e.g., Tremonti et al. 2004, 
iBrinchmann et al.ll2004l) are required to analyse correlations 
over a range in absolute magnitude comparable to the NFGS. 
The NFGS on the other hand, is very well suited as a local 
comparison without corrections and is already sufficient large 
for our purposes. 

5.1.1. Observed correlations with My 

We will proceed by first quantifying the correlations with 
My of indicators which might be responsible for the blue 
CMR, most importantly dust-reddening, gas-phase metallic- 
itym and stellar age. Then we investigate whether a self- 
consistent solution can be found. 

Figure [8] shows the U-V color versus absolute V-band 
magnitude of NFGS galaxies, revealing the well-known color 
bimodality of field galaxies (e.g., Tully et al. 1982, Baldry 
et al. 2004). A blue and a red color-magnitude relation of 
roughly similar U-V slope are visible. A linear fit the blue 
CMR results in 

U-V = -0.078(±0.014)Mv- 1.15(±0.05) (5) 

and can be traced from My = -16 to My = -22. The to- 
tal systematic color variation in this range is « 0.5 mag. We 
define blue sequence galaxies as those having colors within 
A(f/-y) < 0.5 from the best-fit Hnear relation (106 out of 
196) galaxies. From hereon we will analyze the spectral 
properties only of the blue sequence galaxies. To ensure 



reliable measurements of Ha/Hf3 line ratios we limit our- 
selve s to the 91 sources with raw Wnp > 3 (see iJansen et al.1 
|2000|) . Seventy-five procent of the remaining sample have raw 
Whp > 5 a and Whu > 21A. This quality cut removes a few 
more luminous galaxies than fainter galaxies with low Wna, 
but not to a degree that significantly biases the derived rela- 
tions or analysis. We note that the main reason for the absence 
of low Whu at the highest luminosities is the color cut, which 
removed the passively evolving early-type galaxies. 

Figure |9ja) shows the dust reddening towards H II regions 
derived from the ratio of integrated Ha and H/3 fluxes and 
assuming the intr insic ratio 2.85 of case B recombination 
(IOsterbrocklll989l ). The fluxes are corrected for Balmer ab- 
sorption and Galactic reddening jJansen et al. 2000). Clearly, 
the color excess towards H II regions correlates strongly with 
My such galaxies more luminous in the V-band suffer more 
attenuation. The linear fit is 

E{B-V)Hn = -0.060(±0.007)My- 1.06(±0.016) (6) 

consistent the trend seen f or star forming galaxies in the SDSS 
(IBrinchmann et al.ll2004l) . Individual galaxies are scattered 
0.10 mag RMS perpendicular to the relation. The contribution 
of reddening to the observed broadband U-V colors depends 
on the difference between the gas phase and stellar contin- 
uum reddening, which is a function of the geometry of the 
dust distribution (e.g., Calzetti et al. 2000). The differential 
reddening can be parameterized by 

E(B-Vhar = fEiB-V)Hn (7) 

Studies of nearby and distant blue galaxies indicate a range of 
/ = 0.5 - 1 (e.g., Kennicutt 1992, 1998; Calzetti et al. 1996; 
Bell & Kennicutt 2001, Erb et al. 2003), hence the extinction 
toward H II regions is at least as high as that of the stellar con- 
tinuum. It is conceivable that the value of / appropriate for 
local late type galaxies is different from that of high-redshift 
galaxies, e.g., perhaps the high-redshift blue sequence galax- 
ies are more like the local UV-bright starbursts, which have 
/ ~ 0.5 (Calzetti et al. 2000). However, high-redshift galax- 
ies are also different in certain respects, e.g., they are more ac- 
tively star forming for the same amount of obscuration, even 
taking into account selection effects (Adelberger & Steidel 
2000) and they are likely a more heterogenous sample span- 
ning a broader range in stellar masses and SFHs. Also, the 
values derived for local UV-bright starburst of Calzetti et al. 
(2000) rely on colors and spectra in small apertures. Hence, 
aperture effects should be taken into account when comparing 
to the integrated properties of high-redshift galaxies, which 
could increase /. 

Fig. 9(b) shows the N2 =log([NII]A6584/Ha) index, 
which is a single valued indicator of oxygen abundance. 
The N2 index is also sensitive to ionization parameter 
dKewlev & DopitaH2002l) . but it has been shown to provide 
consistent 0/H abun dances for galaxi es with electron tem- 
perature measurement lYin et al.l (l2006h . The absolute abun- 
dances are uncertain, but for our purposes, we are interested 
in the relative abundances, which should be more robust using 
the same indicator on similar galaxies. N2 shows a strong cor- 
relation where more luminous galaxies are more meta l-rich, 
as ha s been firmly established at low redshift (Tremon ti et al.l 
12004 . Using the calibration of iDenicolo et all (i2002i) 

12 + logiO/H) = 9.12 (±0.05) + 0.73 (±0.10) A^2 (8) 

and 12 + log(0/i / ) = 8.69 as Solar metalHcity 
dAUende Prieto et all l2001h . we find an approximately 
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Fig. 9. — Emission lines properties versus absolute V-band magnitude of the blue sequence galaxies in Fig[8] (a) the dust reddening E{B — V)hii towai'ds 
H II regions derived from the Balmer decrement. The solid line shows a linear fit to the data, (b) The N2 =log([Nll]A6584/Ha) metallicity index. The dotted 
line indicates [N11]/Hq corresponding to Solar metallicity. The dashed line shows a linear fit of [NII]A6584/Ha versus My. (c) The Ha equivalent width Whq 
versus My. The dashed line shows a linear fit to the data without any correction for differential reddening of the stellar continuum relative to H II regions, i.e. 
/ = 1 where E(B—'V)star = f E{B—V)hii- However, any trend of W^q with My depends on the differential dust correction, or the geometry of the dust, (d) W^q 
corrected for a differential reddening of / = 0.64, the value obtained from the best fit to all spectra and photometry (see the Appendix). The dashed line shows 
the linear fit to W/f^ cor- There is no significant remaining trend with of W^a-cor with My. 



linear relation between My and metallicity Z of 

Z/Z0 =-0.22(±0.018)Mv-3.32(±0.25) (9) 

The metallicity ranges from ~ 0.4 Zq at My = -17 to 1 .5 Zq 
at My = -22. The luminosity-metallicity dependence derived 
for the SDSS dTremonti et al.l l2004) is slightly steeper, Ukely 
because of different sample selection and metallicity cali- 
bration; e.g., the N2 indicator used here is known to satu- 
rate above approximately Solar metallicity (" Kewlev & Dopital 
1^02). For our purposes, the true form of the relation is 
not terribly important as long as the trend and amplitude 
are roughly correct, because the effects of metallicity on the 
broadband colors of blue galaxies are small (see the Ap- 
pendix). 

Fig. 9(c) shows the Ha equivalent width Whu after correc- 
tion for Balmer absorption. The Ha equivalent width is an 
indicator of the instantaneous star formation rate per unit op- 
tical luminosity, which is often interpreted as the star forma- 



tion rate per stellar mass (specific star formation rate; sSFR) 
or the ratio of present to past averaged star formation. Fitting 
a linear relation log(WHQ) and My yields 

logiWHa) = 0.032(±0.016)My-i-2.06(±0.22) (10) 
Note that the intrinsic slope of log(WHa) versus My may 
be different from the observed slope in a way that depends 
on the amount of differential extinction towards emission line 
regions and the stellar continua /. This can be understood 
as follows. If both Ha and the underlying stellar continua 
at 6563 A are equally attenuated (/ = 1), then the equiva- 
lent width \og(FHa/F\(,5(,3) remains unchanged. If Ha is 
more attenuated (/ < 1), then the observed equivalent width 
changes by a certain amount (/- 1)E{B -V)hii, where the 
exact amount depends on the level of reddening. Because 
for our sample the reddening depends on absolute magnitude 
(Eq. 6), this means the change in Wna increases with abso- 
lute magnitude, and therefore the slope of Slog(WHa)/SMy 
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changes. The extinction correction that we need to apply to 
the observed slope to obtain the intrinsic slope is thus propor- 
tional to (/- l)5E{B-V)Hnl5Mv. For a lCalzetti et aTf(l2000h 
dust-law the slope reverses sign for / < 0.6 (similar values for 
other dust-laws) 

5.1.2. What is responsible for the U — V versus My relation? 

We now use the observed emission line correlations to ex- 
plicitly calculate the contributions of the stellar continuum, 
the metallicity, and the dust attenuation to the observed U-V 
versus My slope for local blue galaxies. The contributions 
cannot be derived uniquely from the individual spectral indi- 
cators, but we can derive simple relations between them using 
l&ruzual & Chariot ( 2003) stellar population models and solve 
for any unknowns. Details are presented in the Appendix, but 
we will discuss the main assumptions and result here. 

The naive way would be to assume a (presumably) ap- 
propriate star formation history and relate the observables 
through the model parameters. The drawback is that any re- 
liance on star formation history might compromise the robust- 
ness of the outcome. For example, even relations that are con- 
sidered fairly robust, such as the relation between the Whc 
and the specific star formation rate still depend substantially 
on assumed SFH. The variation in derived sSFR between ex- 
ponentially declining SFHs with timescales t = 1-30 Gyr is 
approximately ^ 0.4 dex (a factor of 2.5) over the range of 
Whu spanned by our data (10-lOOA). 

Instead we will express our observables in terms of the dust- 
free U-V color of the stellar continuum at fixed metallicity 
iU - V)siar- The main advantage is that our observables, the 
correlations of observed U-V, Whu, E{B-V)Hn, and Z with 
My are easily expressed in (U - V)star with very little depen- 
dence on assumed SFH. The price is that we cannot interpret 
any trend of (U -V)srar directly in terms of stellar "age". As 
shown in the Appendix, we form a linear system and solve 
for the unkown S{U-V)star/ SMy and the ratio of gas-phase to 
continuum absorption /. 

To summarize, the observed correlations are: 



Sy(U-V)ob., =- 
6yE(B-V)Hn = 

6yZ = 



-0.078(±0.014) 
-0.0060(±0.007) 
0.032(±0.016) 
-0.22(±0.0018) 



where Sy denotes S/SMy, the slope with respect 
frame My. 
The unique solution is: 

/= 0.64(±0.14) 

Sy(U - V)„ar = -0.005(±0.016) 



The contributions to the observed U- 
the blue CMR break down as: 



V versus Mi 



Sv{U-V)dus, = -0.073 ±0.016 
SviU-V)„ar = -0.004±0.015 
6yiU-V)z = -0.0014 ±0.003 



(11) 
(12) 
(13) 
(14) 

to rest- 



(15) 
(16) 
(17) 
slope of 

(18) 
(19) 
(20) 



We note that'Brinch mann et alj 120041) find that ignoring the metallicity- 
dependence of the Case B Ha/H/3 ratio would lead to an overestimate of the 
dust attenuation for the most metal rich star forming galaxies, whereas a fixed 
conversion factor from Ho to SFR would underestimate the star formation 
rates. By coincidence, these effects nearly cancel when calculating the dust- 
corrected Ho equivalent width, leaving the results here almost unaffected. 



The result is that only a trend of increasing dust-reddening 
with optical luminosity explains the spectra and broadband 
U-V colors self-consistently, with minor contributions to 
the slope from the dust-free U-V stellar continuum (at fixed 
metallicity) and metallicity. The solution is unique and sys- 
tematic variations are smaller than the random uncertainties 
for all metallicities Z/Zq = 0.2-2.5 and SFHs with e-folding 
timescales r > 1 Gyr. 

The best-fit differential absorption / = 0.64(±0. 15) is inter- 
mediate between those published in literature (e.g., Kennicutt 
1998; Calzetti et al. 1996; Bell & Kennicutt 2001), and it 
implies //// regions suffer up to 0.5 mag more V-band a ttenu- 
ation than the stellar continua for a lCalzetti et al.l (120001) dust 
law. 

There is virtually no variation of the dust-free U-V stellar 
continuum iU-V)s,ars over 5 magnitudes of My. This can be 
seen in Fig. 9(d) which shows the absorption corrected equiv- 
alent width of Ha, which scales directly as SylogWna.cor = 
0.^36y(U-V)stars and shows no remaining trend with absolute 
magnitude. The small variation in the dust-corrected Whq 
would seem to imply a constant SFR/M^ with luminosity, 
in contrast to what is derived for galaxies in th e SDSS from 
optical emission lines ('Brinchmann et al."'2004') and GALEX 
near-UV photometry (Salim et al. 2007). However, we note 
that the dust-corrected Whq measures the line flux (SFR) per 
unit continuum luminosity, which is not the same as SFR/M*, 
because it strongly depends on the assumed star formation 
history (see Fig. Al). Our modeling therefore does not place 
direct constraints on the specific star formation rate SFR/M^,, 
but instead it attempts to determine which trend contributes 
most to the emergent U-V color versus My relation. 

While there is no significant trend of \og(WHa,cor) with 
luminosity, there is substantial scatter at fixed My. The 
scatter is <j(}og(WHa,cor)) = 0.20. Using the model relation 
<T{log(WHa) = 0.83cr(t/- V)j,ar to convcrt fZ-V stellar contin- 
uum colors (see the Appendix), we find a(U -V)srar = 0.24. 
Hence, variations in sSFR are probably the dominant contri- 
bution to the observed scatter (j{U -V)ohs = 0.25. However, 
we cannot exclude that random variations in / also contribute 
to the scatter 

To summarize our main assumptions: the star formation 
history of the galaxies averaged over the blue sequence is as- 
sumed to be smooth, which is plausible even if the SFH of in- 
dividual galaxies is more irregular (e.g.. lRudnick et an i2006). 
and the average integrated SFR is expected to be constant or 
declining, as is the case at z < 1 (jBrinchmann et al. 2004). In 
addition, we assume that the ratio of gas-phase to continuum 
absorption f is approximately independent of My, as is com- 
monly done dCalzetti et al.l2000l:lTYemonti et al.l2004 . How- 
ever, we point out that this has not been verified directly, as 
the distribution of dust is difficult to constrain observationally, 
and it should be addressed in future studies. 

Our preliminary conclusion is that the U-V color- 
magnitude relation for blue galaxies is dominated by system- 
atic variations in dust reddening of the stellar continuum, with 
only minor contributions from intrinsic color variations of the 
stellar continuum ("age") or metallicity. On the other hand, 
the U-V color scatter perpendicular to the blue CMR is likely 
dominated by variations in the stellar ages or specific star for- 
mation rates. More supporting data is needed for more defini- 
tive conclusions (e.g., rest-frame UV and IR), but we will use 
these as a working hypothesis for the rest of the paper. 
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5.1.3. The origin of the Blue Sequence 

While the dkect effects of metallicity on the U -V stellar 
continuum are small, it is possible that the trend of increased 
reddening with rest-frame My essentially reflects the mass- 
metallicity relation or luminosity metallicity relation such as 
found in nearby galaxies (Bell & de Jong 2000; Zaritsky, Ken- 
nicutt & Huchra 199 4), at low r edshift (Tremonti et al. 2004), 
and at high redshift (lErb et al.1 '2006). Qualitatively, this can 
be understood as follows. The optical depth of a star forming 
galaxy depends on the gas column density and gas metallicity. 
If both are higher for more luminous galaxies then both will 
increase the dust reddening. Using a very crude model and 
local galaxy correlations of gas-fraction, stellar surface den- 
sity, and metallicity with luminosity. Bell, (.2003.) find good 
agreement between V-band optical depth and the total IR- 
to-FUV ratio for nearby galaxies. Using the same model, we 
would find a slope of 5v{U -V)dust = -0.2 which is about twice 
as steep as the data, but not wildly different, especially con- 
sidering the simplicity of the model. Nevertheless, empirical 
calibrations do not explain the mechanisms that set up these 
correlations in the first place. 

5.2. Constraints from the Evolution of the Blue Sequence 
5.2.1. The Blue Color Magnitude Relation 

Taking the redshift independence of the blue CMR slope 
at face value suggests that the process responsible for the re- 
lation was already in place at z ~ 3, or possibly even z ~ 4 
(iPapovich et al.ll2004l) . A strong contribution of stellar age 
or star formation history to the relation therefore seems less 
likely. The relative ages of galaxies along the sequence would 
evolve considerably over such a range in time and the blue 
CMR slope would steepen measurably with redshift. 

On the other hand, the non-evolving blue CMR slope may 
be consistent with an underlying mass-metallicity relation. 
There is evidence that a mass-metallicity relation was already 
estabhshed at z > 2 (Erb et al. 2006) with a ^ 0.3 dex off- 
set to lower metalliciti es but broadly simila r shape to the re- 
lation at low redshift dTremonti et al.ll2004h . Models where 
less massive galaxies are less efficient at retaining the metals 
produced by star formation, would naturally create a mass- 
metallicity relation at any redshift. Nevertheless, it is unclear 
whether galaxies on the z = 3 blue sequence can be directly 
compared to those at th e z = blue sequence. For example, 
lAdelber^er et alj ( l2005h find that the clustering properties of 
luminous (blue sequence) LBGs at z = 2 - 3 are more compa- 
rable to local early-type galaxies then local late-type galaxies. 
In addition, it is unclear if the evolving densities, dust-to-gas 
ratios, gas fractions, and metallicities of the galaxies would 
result in the same effective reddening relation at all redshifts. 

In the end deep, high signal-to-noise NIR spectra for large 
numbers of high-redshift galaxies are needed to constrain the 
origin of the trends with metallicity, dust-opacity, and star for- 
mation history directly (e.g., Erb et al. 2006a,b; Kriek et al. 
2006b). 

5.2.2. The Scatter around the Blue Sequence 

If the color scatter perpendicular to the blue CMR arises 
largely from age variations in the stellar populations, then the 
redshift independence of the scatter provides clues to the for- 
mation history of blue sequence galaxies, allowing us to con- 
strain their ensemble star formation history. This technique 
has been applied before using the color distribution of pas- 
sively evolving red sequence galaxies (e.g.. Bower et al. 1992, 



van Dokkum et al. 1998), although obviously the exact same 
models cannot be used here as the blue sequence galaxies are 
still vigorously forming stars. 

The most simplistic approach would be to directly trans- 
late the color scatter to a variation in stellar age. Assum- 
ing BC03 stellar population models with constant star forma- 
tion, where the dependence of color on age is (5(2200 - V) = 
0.72 Slog(age) between 7.5 < log{age) < 10.5, we find that 
the observed scatter a(22QQ-V) = 0.5 translates to log(fl^e) = 
0.7, or a factor of 5. The age scatter would not remain con- 
stant with time, however, and would decrease for constant 
SFHs because the relative age differences become smaller 
with increasing mean age of the sample. One way of keep- 
ing the color scatter constant is to continously add new young 
galaxies to the sample. 

To illustrate how adding new young galaxies to a sample 
influences the evolution of the zeropoint and scatter, we sim- 
ulate ensemble star formation histories in a simple model. 
Here, model galaxies start forming at z = 3.5 following three 
different time dependent distribution functions for the forma- 
tion rate Ppit), the rate at which new young galaxies are added 
to the simulation: 1) exponentially declining with an e-folding 
time of 2 Gyr, 2) constant with time, and 3) exponentially ris- 
ing with e-folding time 2 Gyr. The timescale is set to (a some- 
what arbitrary) 2 Gyr, which causes a factor of ~ 10 change 
in formation rate between z = 0.5 - 3. For each formation rate 
Pfit), BC03 models were used to calculate the evolution of 
the colors and luminosities, each time adopting three differ- 
ent star formation history for the model galaxies: a declining 
SFR with e-folding time r = 0.5 Gyr, a declining SFR with 
e-folding time t = 2 Gyr, and constant star formation. We 
added a fixed amount E(B- V) = 0.13 of dust to the colors 
and magnitudes, adopting a ICalzetti et alJ (l2000l) reddening 
law. For simplicity a flat distribution in masses is assumed 
and a magnitude-limited sample is simulated by imposing a 
threshold on the luminosities. More sophisticated modeling 
will be explored in §6. 

The results are shown in Fig. 10. Each column represents 
the evolution for a different distribution Pfitz)- The top row 
shows the evolution of the mean stellar population age, where 
age is defined as the time elapsed since the onset of star for- 
mation, while the second and third row show the evolution of 
the 2200 - V color scatter and 2200 - V intercept color of the 
simulated samples. Models that are wildly inconsistent with 
the data are SFHs with r < 2Gyr, or SFHs with r > 5Gyr 
models with declining rates of newly formed galaxies Pf(^)- 
In general a rising Ppitz) increases the scatter, while slowing 
down the evolution of the zeropoint, while a declining Pfih) 
does the exact opposite. 

While these simple models are merely illustrative, an inter- 
esting result is that even models that are reasonably consistent 
with the data, fail to reproduce the detailed shape of the color 
distribution (see histograms in Fig. 1). The reason is that the 
constant and exponentially declining SFRs always produces a 
tail to very blue colors, while they have a "ridge" to redder 
colors, as shown in the schematic in Fig. 1 1, in contrast to the 
observations. This is easily understood: the continuing forma- 
tion of new young systems adds galaxies with very blue col- 
ors to the sample, while the 2200-y oc log(f) behaviour of the 
declining SFHs causes a pile-up towards redder colors and a 
red ridge. A similar profile can be seen for passively evolving 
red galaxies at low redshift and in clusters (Ivan Dokkum et all 
199 8). Note that allowing galaxies to evolve away from the 
blue sequence does not change this picture in an obvious way. 
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Fig. 10. — Modeling of the zeropoint and scatter of the blue sequence using simple ensemble star formation histories. Model galaxies start forming atz = 3.5 
following three different distribution functions Pf for the formation rate: exponentially decHning with an e-folding time of 2 Gyr (left), constant with time 
(middle), and exponentially rising with e-folding time 2 Gyr (right). The top row shows the scatter in age of the model galaxies (time elapsed since onset of star 
formation). The middle and bottom row show the evolution of biweight color scatter and the median color respectively, assuming three distinct star formation 
histories for the stellar populations: declining SFR and e-folding time r = 0.5 Gyr (dashed), declining SFR with e-folding time r = 2 Gyr (dotted), and constant 
star formation (CSF; solid). The fat gray lines show the best linear fits to the data. The hatched gray line shows the scatter after subtracting the photometric 
uncertainties in quadrature. The increase and then drop of the scatter in the left middle panel indicates the formation of a red sequence. 



In fact, if the oldest and reddest model galaxies are allowed 
to evolve off the blue sequence then the scatter of the blue 
sequence would be reduced. To match the observed scatter, 
more newly formed blue galaxies would need to be added, 
exacerbating the blue wing. 

Clearly, another way of introducing and maintaining scat- 
ter on the blue sequence is to allow complex ("bursty") star 
formation histories, which are expected in models for galaxy 
formation (Somerville et al. 2001; Nagamine et al. 2005) and 
inferred from fossile evidence in lo cal galaxies ( Trag er et al.i 
I2000t iForster Schreiber et all I2003D . Modulation of the in- 
stantaneous star formation rate is a very effective way to 
change the ultraviolet-to-optical color rapidly. It requires 
somewhat more careful modeling as large changes in the 
sSFR are accompanied by large changes in luminosity. Hence 



selection effects can substantially alter the color-distribution 
of a magnitude limited sample. We proceed next by construct- 
ing more complex models to explain the detailed shape of the 
blue CMR color distribution. 

6. CONSTRAINTS OF THE COLOR-MAGNITUDE RELATION ON 
THE STAR FORMATION HISTORIES OF BLUE GALAXIES AT 
Z~3 

We explore how the observed shape of the scatter perpen- 
dicular to the blue color-magnitude relation can be used to 
constrain models for the ensemble star formation histories. 
As shown in Fig. 1, the observed scatter has a peculiar shape, 
with a "ridge" to blue colors, and a wing or tail to red colors. 
We will focus on the redshift range 2.2 < z < 3.2, where most 
galaxies in the observations are on the blue sequence, and 
discuss the implications for a selection of formation models. 
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Our basic assumption is that the color scatter arises mainly 
from age variations among the galaxies. We do not know 
whether this is true for high-redshift galaxies, but this assump- 
tion seems supported by the large variation of Ha equivalent 
width in local blue sequence galaxies (see §5.1) and the wave- 
length dependence of the blue CMR scatter at z 3 (see §4.4). 

6.1. The Models 

We explore a set of 4 idealized models, characterized by 
distinct parameterizations of the star formation rates: 1) 
galaxies start forming at random redshifts, form stars at a con- 
stant rate, and after a certain time cease star formation, 2) 
galaxies start at random redshifts and have an exponentially 
declining star formation rate, 3) galaxies start at z = 10 with 
constant star formation and experience random star bursts in 
which they produce a certain fraction of the total stellar mass, 
4) galaxies start at z = 10 and form stars episodically, cycling 
through periods of star formation followed by a periods of 
quiescence. 

While the chosen set is not exhaustive, the models do span 
the whole color space that z = 2 - 3 galaxies occupy and 
the parameterizations of the SFR either resemble those of- 
ten used in literatur e, e.g., constant or exponentially decUn- 
ing star formation tPapovich. Dickinson. & Fergusoril 120011: 
Forster Schreiber et al.l2004tlLabbe et al.l2005l;IShaplev et alj 
2005h . or are motivated by models for galaxy for- 
mation and inferences from fossile evidence in lo- 
cal galaxies, e.g., bursting or episod ic star formation 
(ISomerville et al.l fe oOU iNagamine et aH l2005t iTrager et alJ 
llOOO; Forster Schrei ber et al.ll2003h . 

Each model is characterized by only 2 or 3 free parame- 
ters with ^^500 or more independent star formation histories 
for each parameter combination. We do not introduce addi- 
tional free parameters, such as distributions in the star forma- 
tion rates, star formation timescales, or burst fractions. Often 
these extra parameters are degenerate with respect to the re- 
sultant color distribution and can not be constrained by our 
data. 

We assume solar metallicv lBruzual & Charlol (l2003b mod- 
els and adopt a Salpeter IMF from 0.1 to 100 Mq to gen- 
erate the colors. For every Monte-Carlo realization we save 
the model parameters, colors, magnitudes, SFRs and masses, 
which are compiled into 4 model libraries containing on the 
order of ^ 10^ unique star formation histories. For each 
parameter combination we construct color distributions that 
reflect the changes in the spe cific star formation rate and 
compare them to the data (see Ivan Dokkum et al] 1 19981 and 
iKauffmann et al. 2003 for similar studies at low redshift). 

A desirable feature of this type of analysis is that fit- 
ting of the observed color distribution resolves, in a sta- 
tistical sense, some of the degeneracies that occur in fit- 
ting the broadband colors of individua l galaxies (see e.g., 
iPapovich. Dickinson. & FergusonI 1200 lb . For example, a 
galaxy observed shortly after a single burst might have simi- 
lar colors as an older galaxy that is forming stars at a constant 
rate. However, bursting and constant star formation histories 
applied to an ensemble of galaxies give rise to very different 
color distributions. 

6.2. Creating Mock Observations 

We need to account for two aspects of the observations be- 
fore comparing the models to the data: the observations are 
magnitude limited, and the observed scatter contains contri- 



butions from uncertainties in the rest-frame colors and from 
variations in dust content. 

A magnitude limit can substantially alter the resulting 
color distribution if individual galaxies in the models evolve 
strongly in luminosity and color. For example, a faint other- 
wise undetected galaxy undergoing a burst of star formation 
will temporarily turn blue and increase in luminosity. The 
galaxy can now enter a magnitude-limited sample and thereby 
change the color distribution. The amplitude of this effect de- 
pends on the relative number of faint galaxies, hence on the 
steepness of the luminosity function (LF). 

To simulate the effect of a magnitude limit, we adopt a rest- 
frame V-band LF with a faint-end slope a = -1.35 derived 
for a / Tj-selected sample of b lu e galaxies U-V < -25 at z = 
2-3 ( Marchesiniet al.ll2007h . ISteidelet al.N l999') find a = 
-1.6 for t he rest-frame UV LF of optically-selected LBGs, 
while Shap lev et al.l (1200 lb derive an even steeper slope for 
the rest-frame optical LF of the same galaxies, based on a 
correlation between observed R magnitudes and R-K, colors 
in their sample. We find no such correlation in our (deeper ) 
data. We will adopt a = -1.35 of Marchesi ni et al.l ([20071), 
noting that the faint-end slope is considerably uncertain, but 
that the results here are not very sensitive to the exact value of 
a. 

We apply the LF in the following way. We convert a 
Schechter luminosity function with M^^^ = 22.7 and a = 
-1.35 dMarchesini et alj[2007h to an approximate mass func- 
tion usin g an average M/Ly =0 .4 appropri ate for ou r blue 
galaxies dForster Schreiber et al J 12004: Shap lev etal] 12005). 
If M/Ly depends systematically on mass, the slope of the de- 
rived mass function would be different. However, the implied 
change from correlations in our sample is smaller than the 
uncertainty in the adopted faint end slope, so we do not in- 
clude this effect. Random values are drawn from the mass 
function, to well below the mass-limit corresponding to our 
magnitude limit and adopted M/Ly, and the properties of in- 
dividual model galaxies (luminosity etc.) are scaled to those 
masses. This approach allows galaxies to "burst into" or "fade 
out" of the sample when the instant (specific) SFRs, luminosi- 
ties, and colors, evolve with time. 

The photometric uncertainties in the rest-frame colors cor- 
relate with rest-frame My, with fainter galaxies having more 
uncertain colors. We approximate this effect by fitting linear 
relations to the color uncertainties as function of My. We then 
apply the luminosity-dependent uncertainties to the model by 
randomly varying the model fluxes within the errors. The me- 
dian error in the rest-frame 2200 -V color is 0.14 mag in the 
HDFS and 0.19 mag for the MS 1054 field. We include scat- 
ter in the dust-reddening adopting a distribution with mean 
E{B-V) = 0.13 and standard deviation 0.06, where E(B-V) 
is required to be positive. These values are appropriate for 
My = -21 blue galaxies in the local universe (see §5.1) and 
we use them as insufficient spectral data exist for direct mea- 
surements of the reddenin g of blue galaxies at high-redshift'"^. 

Using a ICalzetti et alj fcOOQ) reddening law, we added a 

We note that the mean Ay for ^ = 2 - 3 blue galaxies (mainly "U- 
dropouts") in literature is similar, but the reported width of the Ay distri- 
butions are often broader (e.g., Adelberger et al. 2000, Shapley et al. 2001, 
Papovich, Dickinson, & Ferguson 2001). These estimates, however, could 
suffer from the degeneracy between age and dust in the models, which would 
broaden the distribution, or are difficult to compare to because of differences 
in sample selection and survey depth. In addition, the scatter perpendicular 
to the blue sequence should always be smaller than the scatter for an entire 
sample as the systematic trend of E(B-V) with My has been taken out. 
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Fig. 11. — An explanation of the terminology used in describing the 
shape of the color distribution. Different models for the formation of blue 
sequence galaxies give rise to color distributions with a variety of shapes. 
Two model distributions are shown and the main features are labeled with the 
terminology used throughout the text 



scatter to the model colors of a (2200 — V) dust ~ 0.24 and ap- 
plied extinction to the magnitudes as well. Adopting an SMC 
extinction law instead would result in a(22QQ-V)dust ~ 0.18. 
The precise mean level of dust reddening is not so important 
as we will emphasize fitting the shape of the blue sequence 
scatter. The total 2200 -V scatter added to the models is 0.29 
mag for the HDFS and 0.31 mag for the MS 1054 field. Note 
that the observed scatter at 2.2 < z < 3.2 is significantly larger 
at 0.51 ± 0.09 and 0.62 ± 0.09, respectively (see Table 2 and 
Fig. |5]l. We ignore the modest evolution of the blue CMR 
evolution over the redshift interval z = 2.2-3.2 in both fields, 
which would have added w 0.08 mag to the scatter. 

Finally, we draw samples of galaxies from the model li- 
braries matching the observed redshift distribution and im- 
pose magnitude limits of My = -19.5 for the HDFS and 
My = -20.5 for the somewhat shallower MS 1054-03 field. 
The limits correspond to the magnitudes to which the sam- 
ples are complete for all SED types and galaxy colors. 

6.2.1. Fitting models to the data 

Before fitting the mock distributions to the data, we subtract 
from the observed 2200- V distribution the blue CMR slope 
-0.17(My + 21) (see §4.2). Hence, we will focus on the resid- 
ual color distribution perpendicular to the blue CMR intercept 
at My = -21. 

For each model and parameter combination we calculate the 
two-sided Kolmogorov-Smirnov (KS) statistic of the mock 
color and the observed distribution, giving us the probabil- 
ity that both were drawn from the same parent population. 
The KS test is attractive because it is non-parameteric, does 
not require binning of the data, and is relatively insensitive to 
outliers. Howe ver, it is quite sen sitive to the median of the 
distribution (see lPress et al.|[T992l) . which is undesirable. Our 



principal interest is to study the characteristic shape of blue 
CMR scatter, i.e. the blue ridge and red wing as explained in 
Fig. 11, not the precise value of the median, which may be 
sensitive to our assumptions about the IMF and the amount 
and type of dust added. Therefore, we subtract the median of 
both the model and the observed colors before calculating the 
KS statistic, but we will also discuss the results of a direct KS 
test where interesting. 

Depending on number of free parameters, the procedure re- 
sults in a 2 or 3 dimensional map of KS-test probabilities. 
We find the best-fit parameters by multiplying the probabil- 
ity maps of the individual fields and choosing the parameters 
with the highest probability. We note in advance, that due 
to field-to-field variations the combined KS probabilities will 
never exceed the two-sided KS probability caculated directly 
between the HDFS and MS 1054 observations, which is p=0.7 
to a magnitude limit of My = -20.5. 

6.3. Results 
6.3.1. Model 1: Constant Star Formation 

Model 1 assumes galaxies start forming stars at random red- 
shifts, uniformly distributed in time up to a certain maximum 
redshift z,mix- The SFR of an individual galaxy is constant for 
a certain time t^j after which it drops to zero. The model has 
two free parameters: Zmax, sampled in 20 steps of 0.2 between 
3.2 < Znmx < 7, and f,/, sampled logarithmically in 10 steps 
from 0.08 to 3 Gyr. 

Figure [T2l' a.b.c) shows the color-magnitude evolution with 
time for a characteristic star formation history and the distri- 
bution of the the best-fit model and the observations. Note, 
the blue CMR slope -0.17(My + 21) was subtracted from all 
observations in this section. 

The best-fit parameters are Zmax = 4.6 and fj/ = 1 Gyr with 
fairly low KS probability p=0.07. The reason is that the model 
never reproduces the red wing: the skew towards red colors 
near the peak of the blue CMR. Although the model can pro- 
duce red galaxies in abundance, it does so by creating a build- 
up of passively evolving galaxies in a second red peak (see 
Fig 13). This is a natural consequence of the 2200 -V color 
evolution of passively evolving stellar populations: the colors 
evolve rapidly away from the blue CMR and converge at the 
red peak (see Fig 12a,b). Sudden cut-offs in star formation 
thus do not produce a prominent red wing. The model can re- 
produce the small scatter of galaxies colors while on the blue 
sequence. This follows directly from the oc logf color evo- 
lution of constant SFR models: for approximately constant 
star formation histories a large scatter in stellar population age 
causes only a modest scatter in color 

To understand better the behaviour of this model we calcu- 
lated the confidence regions of the free parameters by boot- 
strapping the observed distribution 200 times, each time find- 
ing the best-fit parameters and identifying the KS-contours in 
the original probability map that encompass 68% and 95% 
of the best-fit Monte-Carlo realizations. The results are over- 
plotted as the blue contours in Fig. 13. The best fits can be 
approximated by fj/ « [{tz=2 + tz=3)~hmax]/'2' = tgvg (the dashed 
white line), where the star formation duration tsf equals ap- 
proximately the mean age tavg of the galaxies at z = 2 - 3. The 
best fits therefore correspond to the situation where a minor- 
ity fraction of the galaxies are just stopping star formation 
and moving off the blue CMR. Models with tsf >> tavg fit 
poorly as they have "inverted" profiles; a blue wing consist- 
ing of new young galaxies that keep being added to the sample 
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Fig. 12. — The results of model 1: random formation redshift and constant star formation for a limited time tsf. {a) The track of 2200- V color against age 
is shown for a characteristic model galaxy (z/- = 4, t^f = 1 Gyr) with filled circles at every lOOMyr. The star formation rate is illustrated schematically by the 
gray line. Also shown are the blue CMR intercepts (at My = -21) in the HDFS {star) and MS 1054 field (diamond), (b) The corresponding track of 2200- V 
color against absolute V-band magnitude in steps of 100 Myr (filled circles), (c) The color-magnitude diagram of the best-fit model (black points) and the data 
(filled gray circles) in the redshift range 2.2 < z < 3.2, shown separately for the HDFS (top) and MS 1054 field (bottom). Only galaxies brighter than the absolute 
magnitude cut-off (da.shed line) are included in the fit. (d) The histograms of residual 2200-y colors of the data (gray histograms) and the best-fit model (hatched 
histograms). The best-fit parameters are Zmax = 4.2 and f,y = 0.9 Gyr. 



and a red ridge because of the log(f ) behaviour of CSF models 
(see Fig 12a). Models with tsf « fm-g have a prominent red 
peak, formed by large numbers of galaxies that are past their 
star formation epoch. 

6.3.2. Model!: Exponentially Declining Star Formation 

Model 2 is similar to the model 1, but now each individ- 
ual galaxy has a single exponentially declining star formation 
rate with e-folding timescale r. This model is characterized 
by two parameters: Zma.K and t, where t is sampled logarith- 
mically in 16 steps from 0.08 to 3 Gyr 

Figure [14] shows a characteristic star formation history and 
the fitting results of model 2. The best-fit parameters are 
Zma.x = 4.2 and t,,/ = 0.5 Gyr with low probability p < 0.01 The 
reason for the poor fit is that the color distribution is always 
too broad and generally symmetric around the median. The 
broad profile is caused by the aproximately linear evolution of 
color with time. The confidence regions of the fitting param- 
eters in Fig. 15 show that the model essentially fits the width 
of the scatter The dashed white line is r oc l/a(t/): a larger 
variety in ages at a certain redshift tr (because of a higher 
Zmax) produces more color scatter a(t~), while longer e-folding 
timescales produces slower evolution with time and less scat- 
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Fig. 13. — The confidence regions for the two parameters in model 1, 
obtained from Monte Carlo simulations. The dark and light blue areas mark 
the 68% and 95% confidence regions and the white cross indicates the best- 
fit. Gray histograms show the model 2200— V color distributions for different 
parameter combinations, p = 0.07 indicates the KS probability of the best 
fit. The dashed white line shows the relation t,f Ri t„rg, where targ = [(t-=2 + 
?_3) -f;„,a.v]/2 or approximately the mean age of the galaxies at z = 2 — 3, 
which provides a good approximation to the best-fits. It corresponds to the 
time where a minority of the galaxies are just turning off star formation and 
moving off the blue sequence. The most prominent variation in the model is 
the build up of a second bump of passively evolving galaxies at t^f < tavg ■ 
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Fig. 14. — Same as Figure [T2l for model 2: random formation redshifts and exponentially declining SFR rate with e-folding time r. An example star formation 
history shown in ia,b) has maximum zo = 3.5 and r = 0.5 Gyr The best-fit color distribution is shown in (c,d). The best-fit parameters of model 2 are Zmax = 4.2 
and T = 0.5 Gyr 
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Fig. 15. — The 68% and 95% confidence regions for the two parameters in 
model 2. Colors and symbols same as Fig. 13. The dashed white line shows 
T (X 1 /a(t,), where alt-) is the standard deviation in the age of the sample 
given a maximum formation redshift z. The most prominent variation in the 
model is the width of the scatter The fit is rather poor as model 2 can only 
produce a broad symmetric blue peak. Low formation redshift are favored, 
because of the emerging red tail of galaxies just moving off the blue CMR, 
which resembles the observations. 



ter. The ratio of the two determines the width of the color 
distribution and this is why the best-fit parameters to the ob- 
served scatter are roughly inversely propertional to each other. 
Low formation redshift are favored, because of the emerging 
red tail of galaxies just moving off the blue CMR which re- 
sembles the observations. 

6.3.3. Model 3: Repeated Bursts 

In model 3, all galaxies start forming stars at a fixed z = 10 
in two modes: constant star formation and, superimposed, 
random starbursts. The bursts are distributed uniformly in 
time with frequency n: the average number of bursts per Gyr. 
During a burst, stars form at a constant elevated rate for a fixed 
time thurst = 50 Myr. The stellar mass formed in the burst is 
parameterized as the mass fraction r of the total stellar mass 
formed so far: Mtursi = fM,ot- Thus, there are two free pa- 
rameters in this model: the burst frequency n and the burst 
strength r. We sample n and r logarithmically in 16 steps over 
the range n = 0.1 -4 Gyr"' and r = 0.01 -4. 

Figure[T6]shows a characteristic star formation history. The 
best-fit parameters are n « 0.3 Gyr"' and r « 3 with a prob- 
ability p = 0.16: extremely massive, but relatively infrequent 
bursts. In fact, for all practical purposes the solution is not 
a "repeated burst" model at all, as very few galaxies undergo 
more than 1 burst between 2 < z < 3. It resembles a forma- 
tion model were a galaxy forms 2/3rd of its stellar mass in an 




Fig. 16. — Same as Figure [121 for model 3. Repeated burst models, consisting of constant star formation and superimposed bursts of a certain strength 
^bursi = f^tot and frequency n. All galaxies start forming at a fixed z = 10. An example star formation history shown in (a,h) has n = 1 and r = 1 Gyr"' . The 
dashed line in (a) shows for comparision a constant star forming model (Bruzual & Chariot 2003). Note that a bursting galaxy population reddens more quickly 
than a constant star forming history. The best-fit color distribution is shown in (c,d). The best-fit parameters of model 3 are n = 0.3 Gyr"' and r = 3; massive, 
infrequent bursts. Massive bursts redden the galaxy population and can reproduce a blue ridge and red wing. However, the absolute absolute colors of the model 
are too red as well. 
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Fig. 17. — (left) The 68% and 95% confidence regions for the two parameters in the repeated burst model 3. Colors and symbols same as Fig. 13. The free 
parameters are the burst frequency n and the bursts strength Mi,i„st = r f^iot ■ The dashed white line shows n oc \/r. Infrequent massive bursts are favored over 
frequent smaller bursts as the latter produces a second blue peak of starburst galaxies which is not observed, (right) The same as the left panel, but now the bursts 
are more obscured by an extra magnitude of visual extinction Ay = \. More frequent smaller burst are now allowed. 



initial burst, after which it forms the rest in residual constant 
star formation over an extended time. 

The model reproduces the correct shape of the color dis- 
tribution, i.e., a blue cut-off and red wing. The blue cut-off 



is the result of the high formation redshift and the lack of 
newly formed, very young galaxies. The red wing is caused 
by galaxies in the post-starburst phase, where large numbers 
of red A-stars formed in the burst outshine the blue O- and 
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Fig. 18. — Same as Figure [T2l for model 4. The episodic star formation model marked by periods of star formation followed by periods of quiescence. The free 
parameters are the duty cycle, or the fraction of the total time spent in activity, the length of the active period ta , and the residual fraction of star formation fraction 
in quiescence r^fy = 0.02. The formation redshift is fixed to z = 10. An example star formation history shown in {a,b) has a duty cycle of 0.4 and r,f,- = 0.02. The 
best-fit color distribution is shown in (c,d). The best-fit parameters of model 4 are dc = 0.67, rjyy = 0.02, and ta = 200 Myr Episodic star formation produces 
the key characteristics of the observed distribution: a blue ridge, red wing, and very blue absolute colors. Note that the stellar population, after resuming stai' 
formation, is bluer than before it stopped, which is why the episodic model evolves more slowly in color than CSF, and which allows the ridge of the distribution 
to remain extremely blue despite the z = 10 formation redshift. 
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Fig. 19. — The 68% and 95% confidence regions for two pai'ameters in 
episodic model 4. Colors and symbols same as Fig. 13. We show the duty 
cycle and the duration of star formation at a fixed residual star formation 
fraction r,-jy> = 0.02 The duty cyle is well constrained. Lower values, and 
hence relatively longer periods of quiescence, would form a second red peak 
of passively evolving galaxies. 



B-stars formed in the residual star formation. The main dis- 
agreement is the absolute mean color, which is too red. This 
is partly because of the higher (z = 10) formation redshift, thus 
higher stellar population ages, but also because successive 
bursts substantially redden a stellar population with respect to 
constant star formation. Hence changing to a lower formation 
redshift does not solve the color offset problem. The reason 
for the redder average color is the dominant contribution to 
the light of the post-starburst population. 

The confidence regions in Fig. [17] show that higher fre- 
quency bursts are disallowed because of the formation of a 
second blue peak, which is not seen in the data. The second 
blue peak can be removed by obscuring the starburst phase 
with an extra magnitu de of visual extinction Ai/ = 1 (see e.g., 
ICharlot & Falll l2000l) . which would allow higher burst fre- 
quencies. We obtain similar results if the extra attenuation is 
more than 1 magnitude. The roughly inverse proportional re- 
lation between the best-fit frequency and burst strength shows 
that both parameters contribute comparably to creating the red 
wing; i.e., the model can now produce a red wing through 
more frequent smaller bursts, or through fewer more massive 
bursts. 
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Fig. 20. — The 68% and 95% confidence regions for episodic model 4 as function of tlie residual fraction of star formation. Colors and symbols same as 
Fig. 19. A higher residual fraction of stai' formation in the quiescent state (less contrast between the active and passive phase) allows lower duty cycles. The 
combined 95% constraints are that the duty cycle is more than ~ 40%, and that the contrast in SFR between the high and low phase, i.e. the relative strength of 
the burst, is more than a factor of ~ 5, which corresponds to more than 0.35 dex of scatter in log{SFR) around the mean. 



6.3.4. Model 4: Episodic Star Formation 

In model 4, galaxies start forming stars at a fixed z = 10 and 
subsequently cycle through a periods of active star formation 
and quiescence. During activity they form stars at constant 
rate for some time and during quiescence the star formation 
rate drops for some time tg to some fraction r,/,. of the active 
value. The primary parameter is the "duty cycle", the fraction 
of total time spent in the active state dc = ta/{ta + tg). For sim- 
plicity we implement the episodic model with repeating duty 
cycles of fixed length with random offsets in the phase. The 
three free parameters are the duty cycle dc, the duration of the 
activity t^, and the residual star formation fraction r^jr. We 
sampled all three logarithmically: dc in 24 steps from 0.05 to 
0.95, ta, 16 steps from 0.08 to 2 Gyr, and r^/,- in 10 steps from 
0.01 to 0.8. 

Figure [18] shows a characteristic star formation history. The 
best fit of the model is dc = 0.67, t„ = 0.2 Gyr, and rjy,- = 0.02 
at relatively high probability p = 0.43, indicating the shape of 
the color scatter is well reproduced. Relatively high duty cy- 
cles are favored, with galaxies spending ^ 70% of the time 
in their active state. Episodic models with high duty cycles 
naturally reproduce the characteristic features of the observed 
scatter The narrow blue peak with a blue ridge in this model 
is caused by the high formation redshift and constant star for- 
mation rates. The red wing results from the short inactive 
period, in which galaxies briefly depart from the blue CMR 
without ever moving far away. Incidently, the episodic model 
is the only model that also produces the correct blue absolute 
colors at z = 2-3 interval (2200- V « -1.0, comparable to 
the observations, see Fig. 4). The origin of the blue color is 



a typical aspect of episodic SFHs: when star formation re- 
sumes after the quiescent period, the color of the population 
is bluer than before quiescence (see e.g.. Fig. 18a). The rea- 
son is that the stellar population formed previously has faded 
during the inactive period, thus on resuming activity there are 
fewer intermediate aged stars (10^ years) contributing to the 
UV luminosity, which is then dominated by younger, bluer, 
stars. The color evolution with time is substantially slower 
than constant star formation helping the model galaxies to 
maintain extremely blue colors despite the high z = 10 for- 
mation redshift. The only mismatch to the data is the deficit 
in the models of very red galaxies A(2200-y) > 1.6, galaxies 
that are 1 .6 mag redder than the blue sequence zeropoint. 

The confidence regions of parameters dc and at fixed 
f.ifr = 0.02 (see Fig.[T9b show that the duty cycle is best con- 
strained, with 68% confidence levels at fixed r^/r of dc = 
0.67(±0.10), while the period of the cycle may vary from 
150 Myr to 600 Myr (68%). Fig. 20 shows the confidence re- 
gions as a function of residual star formation rate rj/, . Clearly, 
if the quiescent star formation rates are higher, then the duty 
cycle can be lower because the model galaxies redden more 
slowly between the burst. However, the contrast in star for- 
mation rate between activity and quiescence must be at least 
a factor of 5 (rj/r < 0.2). 

6.4. Discussion 

We have illustrated how the shape of color scatter might 
be used to constrain simple parameterizations for the ensem- 
ble star formation history of galaxies on the blue sequence. 
Before we discuss its implications, we stress that our models 
are on purpose a limited and idealized set and not necessarily 



21 




-24 -22 -20 -18 -24 -22 -20 -18 -24 -22 -20 -IB -24 -22 -20 -18 

My (rest-frame) 



Fig. 21 . — Predicted evolution of the color magnitude distribution for all models in §6. Each column presents one of the models, while each row gives the 
color-magnitude distribution at different redshifts. The models were fit to the observed distribution at z = 2.3 — 3.2, and subsequently evolved forwai'd in time 
using the best-fit parameters. Model 1 (constant star formation with cut-off), produces a strong red peak at lower redshift, but keeps a blue sequence with low 
scatter. Model 2 (exponentially declining SFR) produce a very broad distributions extending to red colors. Model 3 (repeated bursts) produces a blue sequence 
with a red wing, which is offset to redder colors compared to model 1 . In addition, a second blue peak at bright magnitudes is visible, containing galaxies caught 
in the starburst. Note the bursts grow stronger with time, because they are parameterized as Mturst = rM,oi. Model 4 (episodic star formation) produces a color 
magnitude distribution with a blue sequence that has a blue ridge and a red wing, and which evolves very little with time. 



exhaustive. While the models can clearly reproduce some or 
all of the important characteristics of the blue sequence, there 
may be other paramerizations of the SFR or treatment of dust 
that can reproduce the characteristics as well. In addition, if 
the degree of reddening depends on stellar age (e.g.. Char- 
lot & Fall et al. 2000, Shapley et al. 2003) then the picture 
is complicated further as star formation history and dust evo- 
lution can only be constrained simulatenuously. As such our 
results are not definitive or unique. However, it is certainly in- 
structive to see the ensemble color evolution of more complex 
star formation histories parameterized in a simple way. 

6.4.1. Implications for the Star Formation History of Blue 

Galaxies 

Several general conclusion can be drawn from the observa- 
tions of the blue sequence in §3 and §4 and the modeling of 
ensemble star formation histories in §5.2.2 and §6. 

First of all, it is clear from §6.3.1 and §6.3.2 that the sim- 
plest models (e.g., constant SFR and cut-off, and declin- 
ing SFRs), fit the observed color distribution and its evo- 
lution rather poorly. It can be easily seen from the lin- 
ear oc t and power law oc log(0 color evolution of expo- 
nentially declining and constant SFHs, that simple random- 
ized distributions of such templates can never create the cor- 
rect shape of the color distribution, i.e., a blue ridge and a 
red wing. Such simple SFHs are widely used for model- 
ing the broadband SEDs of individual high-z galaxies (e.g. , 
Shaplev et al."200r; 'Papovich. Dickinson. & Fergusonl l2001l; 
van Dokkum et al. 2004; Forster Schreiber et al. 2004), hence 
their derived properties should be interpreted with caution. 
Nevertheless, it is remarkable that a simplistic CSF model can 
roughly reproduce some of the observed trends with redshift. 



such as the small evolution of the scatter with redshift and 
the approximate linear evolution of the intercept over z = 1 - 3 
(see Fig. 22). This may indicate that blue sequence galaxies 
are actively forming stars at nearly constant rate most of the 
time, and that their SFHs are rather similar over long periods 
> 1 Gyr 

iPapovich. Dickinson. & FergusonI (1200 lb already argued 
that constant or declining SFH are probably wrong for most 
high-redshift galaxies, as they lead to inferred ages 100- 
300 Myr for z = 2-3.5 LBGs in the HDFN with relatively 
small scatter. Such an age distribution, much smaller than 
the cosmic time span of 2 < z < 3.5 (~1 Gyr) is implausible 
unless the galaxies are evolving in sync (impossible) or the 
assumed SFH and hence the implied age is wrong (likely). 
Episodic or bursty star formation histories provide an alterna- 
tive by rejuvenating the appearance of the galaxies, possibly 
allowing much higher ages and inferred formation redshifts. 

The repeated burst model with underlying constant star for- 
mation (model 3) reproduces two key features of the color 
distribution: relatively small scatter, caused by the constant 
star formation, and a red wing, caused by the red colors of 
the post-starburst stellar populations - off setting the model 
galaxies to the red of the blue CMR . The bursts have to be 
massive with mass fraction formed in the burst > 40% be- 
cause the fading starburst needs to outshine the on-going star 
formation. This effect can be clearly seen in Fig. 16, where 
galaxies burst briefly to high luminosities and ultrablue colors, 
then fade to redder colors, and ultimately return to the blue se- 
quence driven by the on-going star formation. The predicted 
evolution of the model is rather different from the observa- 
tions however, with a strongly increasing scatter with time 
(see Fig. 21 and Fig. 22). This is caused by our choice of pa- 
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Fig. 22. — Prediction evolution of the scatter and intercept color of 
the blue CMR in the models of §6: model 1, constant star formation with 
cut-off (blue solid line); model 2, exponentially declining star formation 
(green dotted line); model 3, repeated bursts, (red dashed line); and model 4, 
episodic star formation (black dash— dotted line). The fat gray lines show the 
best linear fits to the observed evolution with redshift. The hatched gray line 
shows the scatter after subtracting the photometric uncertainties in quadra- 
ture. 



rameterization for the burst strength, which was expressed as 
a fraction of the total stellar mass formed, leading to burst that 
are increasingly powerful. The strength of the burst is rather 
problematic, even at early times. At z = 3 the model predicts 
20-80 fold increases in SFR during a burst. This would pro- 
duce blue galaxies with SFRs exceeding 1000 Moyr"', which 
are not observed. That being said, for a frequency < 1 Gyr~' 
and burst duration of 50 Myr, we would only expect approx- 
imately 5 such sources in our fields, hence their absence is 
only marginally significant. In addition, star formation and 
obscuration are known to correlate at low and high redshift 
(e.g., Adelberger & Steidel 2000), so if the starburst episode is 
highly obscured the galaxies might briefly resemble extreme 
heavily obscured starbursts, such as those detected in the sub- 
mm (e.g.. Chapman et al. 2000). 

The episodic star formation model (model 4), here mod- 
eled as a simple step function modulation of the SFR, explains 
the key features of the blue sequence surprisingly well. The 
model reproduces the correct shape of the color scatter, with a 
narrow blue ridge, because of constant star formation and high 
formation redshift, and a red wing, resulting from the short 
inactive period in which galaxies briefly depart from the blue 
CMR. The predicted evolution in zeropoint and scatter (see 
Fig. 22) is roughly consistent with the data. The most impor- 
tant parameter in this model is the duty cycle, or the fraction of 
time spent in the "active" state, which our observations con- 
strain to dc = 0.67 ±0.1 if the contrast in star formation rate 
is a factor of 50. The combined constraint on duty cycle and 
variation of star formation rate is that the observations require 



at least a 40% duty cycle and at least a factor of 5 variation 
of the star formation rate (at 95% confidence, see Fig. 20). A 
factor of 5 in SFR in this binary active-passive model, corre- 
sponds to 0.35 dex of scatter in \og(SFR) around the mean. 
Interestingly, the color evolution for this model is slower than 
that of constant star formation, making the average galaxy sig- 
nificantly bluer for the same age than a galaxy with a constant 
star formation (see Fig. 18a). If blue sequence galaxies at 
Z = 3 really have such episodic SFHs, many of them could 
have started forming at much higher redshifts than would be 
inferred from their average stellar population age, suggesting 
a blue CMR might exist at significantly earlier times. 

6.4.2. Implications for Red Galaxies 

Shapley et al. (2005) have suggested that z > 2 galax- 
ies may pass in and out of UV-selected and near-IR color- 
selected samples as they evolve from phases of active star for- 
mation, when they are blue, to quiescence, when they are red, 
and back again. An interesting question is to see whether an 
episodic model can explain the tail to very red colors of galax- 
ies with A(2200-y) > 1.6 mag. 

In our z > 2 observations red galaxies are almost exclu- 
sively C/V-faint Distant Red Galaxies, which can be selected 
by J,-K, > 2.3 (Franx et al. 2003; see Fig. 1). Half of the lu- 
minous Mv < -20.5 DRGs at z = 2-3 have A(2200-y) > 1 .6 
mag. By number they make up about 15% of all luminous 
My < -20.5 galaxies. We find that our best-fit episodic model 
does not produce substantial numbers of such very red galax- 
ies (only < 3% of My < -20.5), due to the high duty cycle 
and the brief time spend in quiescence. Naturally, an episodic 
model can produce more red galaxies if a minority of the 
blue galaxies are allowed to be passive longer For the red- 
dest galaxies, however, the time needed to passively evolve 
to such colors would approach the Hubble time, so it is ar- 
guable whether such formation histories can be considered 
"episodic". 

Therefore, the model suggests that roughly half of the 
DRGs to Mv < -20.5 could participate in the episodic SFH, 
and thus in principle be blue sequence galaxies in a quiescent 
state. The rest of the DRGs may have turned off star forma- 
tion altogether, as suggested by analysis of the IRAC imaging 
in the HDFS, or dust obscured star formation may play a role 
in causing the red colors (e.g., Labbe et al. 2005). Recent 
studies of DRGs indicate that half of them are detected with 
Spitzer/MIPS and are forming stars at a fairly high rate while 
heavily obscured (van Dokkum et al. 2004, Papovich et al. 
2006, Knudsen et al. 2006, Webb et al. 2006). Hence dif- 
ferent mechanisms may be at work to cause the red colors. 
Interestingly, the high detection rate with MIPS of DRGs also 
suggests high duty cycles, similar to what we infer for blue 
galaxies. 

Currently, we fitted only the color distribution as our sam- 
ple is too small to split up in luminosity. Models where red 
galaxies are passive phases of blue sequence galaxies pre- 
dict that the characteristic V-band luminosity of red galax- 
ies is 1 - 2 mag fainter than that of blue galaxies. There is 
no clear evidence in larger somewhat sh allower multiwave- 
length surveys that this is the case (e.g., Ivan Dokkum et al.1 
2006; Marchesini et al. 2007). With future larger and deeper 
samples, we can attempt to model the color and luminosity 
evolution more consistently. 

6.4.3. Summary 
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To conclude, we briefly summarize the main characteris- 
tics of the blue sequence color distribution, and what type of 
models might explain them. 

Relatively small scatter, constant with redshift is naturally 
reproduced by model galaxies that are forming stars at a 
(nearly) constant rate on long timescales, suggesting that 
galaxies spend a long time on the blue sequence in active for- 
mation, or equivalently, have high duty cycles. Any model 
with declining star formation rates and t < 2Gyr fits the data 
badly. 

A Red wing is hard to produce with a smooth, slowly vary- 
ing SFHs, as their colors evolve with time as cx f or oc log(f ); 
this always produces a red ridge (see Fig. 11). Pure dusty 
origin of the red wing is unlikely, as the predicted scatter in 
far-ultraviolet colors is larger than observed (see Fig. 6). The 
red wing can more easily be created by a temporary drop in 
star formation rate, suggesting that galaxies form with non- 
monotonic SFRs. Simple episodic models with a binary high 
and low level of star formation can reproduce the amount scat- 
ter, its shape, and its evolution. The best fit-parameters are a 
high duty cycle, where model galaxies spend 70% of the time 
in the active star forming state, and a variation in SFR between 
the active and passive phase of a factor of 50. Duty cycles less 
than 40%, and a contrast in SFR less than a factor of 5, corre- 
sponding to 0.35 dex of scatter in \og{SFR) around the mean, 
are excluded at the 95% confidence level. We caution that 
our working hypothesis that the scatter of the z = 2 - 3 blue 
sequence is mainly attributable to variations in specific SFR 
must be confirmed by deep NIR spectroscopy (e.g., Kriek et 
al. 2006a,b). 

The blue ridge puts limits on the production rate of newly 
formed galaxies with very young stellar populations, or limits 
on the frequency of massive star bursts. The exact limits are 
model dependent, but the existence of a blue ridge at z = 3 sug- 
gests that the majority of the blue sequence galaxies arrived 
on the sequence long before the epoch of observation, i.e. for- 
mation redshifts subtantially higher than z = 3. Powerful star 
bursts (leading to a 20-fold or more increase in SFR) must 
occur at low frequency < < 1 Gyr"' or be heavily obscured. 
Otherwise the combination of luminosity increase and a steep 
luminosity function would create a second starburst peak in 
the color distribution, blueward of the blue CMR (see Fig. 17 
and Fig. 21), which is not observed. 

Blue absolute colors, evolution with redshift can be ex- 
plained by episodic star formation, rising star formation rates 
before z = 3, or a different "top-heavy" IMF, which produces 
more blue 0/B stars per unit stellar mass formed. Such con- 
stant/rising/episodic SFHs cannot continue to the present day, 
as they predict very modest evolution of the mean color at 
z < 1, resulting in too blue colors for the blue sequence at 
z = 0. The well-established decline of the global SFR be- 
tween z = 1 and z = (e.g., Lilly et al. 1996, Cowie et al. 
1999, Glazebrook et al. 1999) does not necessarily contradict 
a more modest decline of the average SFR of blue sequence 
galaxies, as much of the global decline is caused by galaxies 
stopping star formation and evolving off the blue sequence al- 
together (e.g.. Heavens et al. 2004). Nevertheless, the evolu- 
tion of the average blue CMR color at fixed luminosity would 
seem to require that blue sequence galaxies are also gradually 
slowing down star formation z < 1 . Possibly, the evolution of 
the mean dust reddening contributes as well to the redder col- 
ors, as it may follow the trend of metallicity which is observed 
to increase from z = 3 to z = (lErb et al..,2006.) . 



7. THE ONSET OF THE RED GALAXIES 

The color distribution in the FIRES fields evolves strongly 
over z = 1-3, especially at the bright end My < -20.5 (see 
e.g.. Fig. 1). At z 3 most galaxies are on the blue sequence 
and there is no evidence for a well-populated red peak. The 
onset of a red peak is tentatively observed in the histograms at 
z= 1.5-2 in the field of MS 1054, whereas it is clearly visible 
at z 1 . Prominent red sequences in the field are also seen in 
photometrically selected samples up to z = 1 (e.g.. Bell et al. 
2004b, Kodama et al. 2004). 

We use a simple color criterion to define and quantify the 
redshift evolution of the red galaxy fraction. In the litera- 
ture, color-limited fractions are usually defined rel ative to the 
red color-magnitude relation (e.g.. Butc her & Oemler 1984). 
However, such a definition is unusable here as the red se- 
quence is virtually absent in our sample at z ~ 3. Hence, we 
define red galaxies relative to the blue sequence as any galaxy 
that is A(2200-y) > 1.6 mag redder than the best fit blue 
color-magnitude relation (see §4.3). The threshold is about 
3 times the rms scatter of the blue CMR and corresponds to 
approximately A(U-V) > 0.7. Note that there is no evidence 
from the data that the scatter depends strongly on redshift. 

Figure |23] shows the evolution of the red galaxy fraction 
by number (Nred/N,„,) and by luminosity density jv.red/ jv.ioi- 
We calculated the fractions in both fields to a fixed rest-frame 
magnitude limit My = -20.5, to which we are complete for all 
SED types to z ~ 3. The luminosity densities were computed 
by adding the luminosities of the galaxies above the magni- 
tude limit, and dividing by the cosmic volume of the redshift 
bin. Uncertainties are obtained with bootstrapping. 

In both fields, we find a sharp increase in the fraction of 
red galaxies from z ~ 3 to z ^ 0.5 and the fractional V-band 
luminosity density. A linear fit in time fits the data points 
reasonably well and we find a red fraction fred that increases 
as 

/„^=-0.07(±0.09) + 0.08(±0.017)f (21) 

with t the time of the universe in Gyr. 

Field-to-field variations are substantial. In the z = 0.7- 1.4 
bin much of the difference can be explained by the contribu- 
tion of the cluster at z = 0.83 in the field of MS 1054-03, which 
contains many bright early-type galaxies dvan Dokkum et all 
I2OOOI) . Eliminating those from the sample - by removing the 
volume between 0.81 < z < 0.85 with mostly spectroscop- 
ically confirmed cluster members - shows they contribute 
about half of the red galaxies. This difference is compara- 
ble to the field-to-field variation at other redshifts (a factor of 
^2). 

The strong evolution of the red galaxy fraction at high 
redshift contrasts to the modest evolution at z < 1. Bell et 
al. (2004b) find a constant luminosity density in photometri- 
cally selected red galaxies in the range 0.2 < z < 1. While 
their selection criteria were different (more stringent), it does 
suggests that much of the evolution in the red galaxy popula- 
tion occurs between z = 1 and z = 3. Recent studies i n other 
fields have suggested the same jGiallongo et al.ll2005h . How 
the increase in the red galaxy fraction relates to formation of 
the passively evolving early-type galaxies is not well known. 
In the low redshift universe, ^ 70% of the red peak galaxies 
from the SDSS are morphologically early-type (Strateva et al. 
2001; Hog g et al . 2002). Data from the COMBO-17 and 
GEMS ( Rix et al.ll2004 , )_ survey indicate this picture is largely 
unchanged to z 0.7 dBell et alj|2b04al) . but the majority of 
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Fig. 23. — Evolution of tlie fraction of red galaxies as a function of redshift. in the HDFS (filled circles) and MS1054 (diamonds). In (a) we show the 
number fraction of galaxies that are A(2200-y) > 1.6 mag redder than the blue color-magnitude relation. This color cut isolates galaxies that belong to the blue 
sequence and corresponds to ~ 3 X the scatter on the blue CMR (see Fig. 5). The dashed line is a linear fit to the number fraction as a function of cosmic time. 
(b) The relative fraction of the total luminosity density in red galaxies. There is a significant increase in the red galaxy fraction between z = 3 and z = 1. We note 
that the field of MS 1054 contains a massive cluster at z = 0.83. The squares in (a) and (b) show the sample after removing the volume between 0.81 <z < 0.85, 
which contains mostly spectroscopically confirmed cluster members. 



the red galaxies at 1 < z < 3 show signs of dust-reddened 
star formation (e.g., Yan & Thompson 2003; Moustakas et al . 
2004, [ van Dokkum e t al.ll2004l: lF5rster Schreiber et al.ll2004t 
iPapovich et al.<,200Q) . 

We defer such interpretations for our sample, until we bet- 
ter understand the nature of the red galaxies at z > 1. The 
foremost questions are which fraction of all red galaxies are 
truly early-types with passively evolving stellar populations, 
and what other factors, such as reddening by dust play a 
role. Rest-frame optical morphologies from space and deep 
mid-infrared imaging with Spitzer/IRAC are providing some 
of the answers ( e.g., Labbe et al. 2005: Pa povich et alJl2006l : 
IZrrm et an2006l) . although large samples are very hard to ob- 
tain. 

The next step in this kind of analysis would be to estab- 
lish at what redshift the narrow red sequence establishes itself. 
Unfortunately, our photometric redshifts are too uncertain to 
allow a determination rest-frame colors with an accuracy bet- 
ter than 0.04 mag, typical of the sca tter in the local red color- 
magnitude relation (Bow er. Lucev. & Ellis 199^. Hence NIR 
spectroscopy is needed to establish the onset of the red color- 
magnitude relation. 

8. SUMMARY AND CONCLUSIONS 

We used deep near-infrared VLT/ISAAC imaging to study 
the rest-frame color-magnitude distribution of infrared se- 
lected galaxies in the redshift range 1 < z < 3 and compared 
their properties to galaxies from the NFGS at z = 0. Con- 
trary to the situation at z = 0, where the color distribution has 
a prominent red and blue sequence of galaxies, we find no 
evidence for a well-defined red sequence at z ^ 3. We did 
find a well-defined blue sequence of star-forming galaxies at 
all redshifts. These galaxies populate a blue color-magnitude 
relation, such that more luminous galaxies in the rest-frame 
V-band have mildly redder ultraviolet- to-optical colors. 

The slope of the blue CMR in our fields is constant up to z ^ 



3, with S(U-V)/SMv = -0.09 ±0.01 or S(2200-V)/5Mv = 
-0.17 ±0.02; identical to the slope of blue, late-type galax- 
ies in the local universe. Using the observed correlations of 
spectral indicators with absolute V-band magnitude of nearby 
galaxies from the NFGS, we explicitly calculated the contri- 
bution of the stellar continuum, dust reddening, and metal- 
licity to the blue color-magnitude slope. Assuming the ratio 
of stellar continuum reddening to gas-phase reddening does 
not depend on absolute V-band magnitude, we find that the 
local blue CMR can be explained almost exclusively by the 
trend of increasing dust reddening with increasing absolute 
V-band luminosity, with only minor contributions from the 
stellar continuum and metallicity. 

While the slope of the blue CMR remains constant up to 
high redshift, the zeropoint at a given magnitude reddens 
gradually from z ~ 3 to z = 0. The evolution over this red- 
shift range is A{U-V) 0.75 or A(2200- V) « 1.4 and the 
reddening is approximately linear with redshift. Most of the 
evolution can be explained by aging of the stars, but the pre- 
cise form of the average star formation history is uncertain. 
The very simplest star formation histories, constant and ex- 
ponentially declining, do not reproduce the details of the ob- 
served color evolution very well. Such formation histories 
assume that the galaxies remain on the ridge of the blue CMR 
throughout their life, which may very well be wrong. Another 
possibility is that evolution in the dust content contributes to 
the evolution colors. However, age-dust degeneracies make it 
difficult to seperate the contribution of age and dust. 

A key feature of the blue CMR relation is that the color 
scatter is markedly asymmetric, with a blue ridge and a wing 
towards red colors. Assuming the scatter is caused by varia- 
tions in the specific star formation rate, we have constructed 
models to explore the constraints that the width and shape of 
the scatter place on the star formation history of blue sequence 
galaxies at z = 2-3. We find that models where galaxies on the 
blue sequence only have simple constant or declining star for- 
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mation rates do not fit well. However, models with episodic 
star formation provide a good fit to the width and shape of the 
scatter. In this model, galaxies form stars with a SFH char- 
acterized by alternating periods of active star formation and 
quiescence. The combined constraints on the best-fit parame- 
ters of the episodic model are that the duty cycle, the fraction 
of time spent in the active state, is more than 40%, and that the 
contrast in SFR between the high and low phase, i.e. the rela- 
tive strength of the burst, is more than a factor of 5, which cor- 
responds to more than 0.35 dex of scatter in log(SFR) around 
the mean. The z = 2-3 model galaxies spend on average more 
than ~ half of the time actively forming stars, while for the re- 
mainder they form stars at a reduced rate and briefly evolve to 
redder colors, forming the red wing. The period of the duty- 
cycle is poorly constrained by the data as it correlates with 
the amount of residual star formation during the quiescent pe- 
riods; the data allow a range of 150 Myr to 600 Gyr. 

Interestingly, the absolute colors of episodic models are 
bluer than those of constant star forming models of the same 
age and evolve slower with redshift as the colors are rejuve- 
nated periodically. Hence episodic models allow higher ages 
and formation redshifts for the same color and do not lead to 
the somewhat enigmatic results from previous studies, where 
the best-fit ages of blue z 2-3.5 galaxies are always much 
smaller than the time interval spanned by the observed red- 
shift range (e.g., Papovich, Dickinson, & Ferguson 2001). 

The episodic models analyzed here do not produce substan- 
tial numbers of very red galaxies, A(2200-y) > 1.6 mag or 
A{U -V) = 0.7 mag redder than the blue CMR. In our data 
these red galaxies are almost exclusively Distant Red Galax- 
ies, selected by observed J, -Kg > 2.3 (Franx et al. 2003). 
Shapley et al. (2005) have suggested that z> 2 galaxies may 
pass in and out of UV-selected and near-IR color-selected 
samples as they evolve from phases of active star formation 
and back again. The results here indicate that half of the 
DRGs are probably too red to be blue galaxies in the quies- 
cent phase of the duty cycle. These galaxies may have turned 
off star formation altogether, as suggested by analysis of the 
IRAC imaging in the HDFS (e.g., Labbe et al. 2005). Alter- 
natively, the reddest galaxies are very dusty and star forming 
at high rates, as indicated by NIR spectra (van Dokkum et 
al. 2004) and IR/Sub-mm observations (Papovich et al. 2006, 
Knudsen et al. 2006, Webb et al. 2006). 

While the modeling of complex star formation histories is 
instructive, we caution that the choice of models and SFR pa- 
rameterizations presented here remain a limited and idealized 
set. As such we can not be certain that the results are unique. 



For example, we have not considered effects that plausibly 
contribute to the evolution of the colo r magnitude distribu tion, 
such as merging of galaxies (see e.g.. lBaldrv et al.ll2004l) . 

Finally, we find a sharp increase (a factor of 6) between red- 
shift z = 2.7 and z = 0.5 in the relative number and the relative 
rest-frame V luminosity density of luminous (My < -20.5) 
red galaxies, where we define red galaxies as redder than 
A(2200-y) > 1.6 relative to the blue CMR. Studies at red- 
shifts z < 1 imply little evolution in the red galaxy population, 
suggesting that the bulk of the evolution takes place between 
z ~ 3 and z 1. There is substantial variation (x2) between 
the red galaxy fraction in the HDFS and MS 1054-03 fields in 
all redshift bins. Obviously, our fields are still very small and 
the variations between the fields suggest that the uncertain- 
ties in the red galaxy fraction are dominated by large scale 
structure. Future sample with very deep NIR imaging will 
provide better estimates. While we are likely viewing the on- 
set of the red passively evolving galaxies that will form the 
well-known red sequence at lower redshifts, additional spec- 
tra and imaging at longer wavelengths is needed to identify 
which of the red galaxies are passive or star forming. Such 
studies are even becoming possible for galaxies too faint for 
spectroscopy, thanks to the infrared IRAC and MIPS instru- 
ments aboard the Spitzer Space Telescope (see e.g. Labbe et 
al. 2005, Papovich et al. 2006). 

On the other hand, the evolution of the slope, zeropoint, and 
scatter of the blue color magnitude relation is similar between 
the fields, indicating that field-to-field variations do not play a 
large role there. Here the challenge is to confirm directly the 
cause of the relation between color and magnitude, and the 
cause of the blue ridge and skew to red colors of the scatter 
around the relation. Very high signal-to-noise spectroscopy in 
the near infrared will be required to measure the balmer emis- 
sion lines in order to estimate ages and reddening. Further- 
more, these studies need to be performed at higher redshifts. 
The advent of multi-object NIR spectrographs on 8-lOm class 
telescopes may make such studies feasible in the near future. 
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APPENDIX 



THE ORIGIN OF THE BLUE COLOR-MAGNITUDE RELATION IN THE LOCAL UNIVERSE 

We use the observed correlations of spectral indicators with absolute V-band magnitude of nearby galaxies to explicitly 
calculate the contributions of age, dust, and metallicity to the blue U -V color versus V-band magnitude slope. Writing 
Syx to denote the slope of quantity x versus absolute V-band magnitude, we want to know Sv{U -V)age, <5v(t/- V)du.«, and 
dv(U-V)z- The quantities cannot be derived uniquely from the spectral indicators or broad band colors separately. However, the 
combination of the broadb and color slope and the spectral observables allows for a unique solution with few assumptions. Using 
iBruzual & Chariot (l2003h stellar population models, we derive simple relations between the spectra and colors. We express 
the relations in terms of two free parameters: the unknown dust-free U -V slope of the stellar continuum at fixed metallicity 
Sv(U -V)s,ar, and the ratio of gas-phase to continuum absorption /. Then we form an appropriate linear system and solve for 

iU-V)„armdf. 
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(u-v),t„ (U-V),t„ 

Fig. 1. — (a,h,c) Relations between vai ious iBruzual & Charloll 120031) stellar population parameters and the dust-free equivalent width of Ha (Whol)- The 
tracks are 4 star formation histories, from top to bottom: constant star formation, and exponentially declining star formation with e-folding times t = 5,2, 1 Gyr 
{a) Whci versus the age, the elapsed time since the onset of star formation, (h) Wfja versus specific star formation rate (sSFR) or SFR/M, . (c) Vfua versus the 
dust-free color of the stellar continuum. The relation to age, sSFR, and color is increasingly better approximated by a linear fit, with decreasing dependence on 
assumed SFH. (d) The difference mU — V color beteen metallicities Z = 0.04 and Z = 0.05, where Z = 0.02 is Solar Two star formation histories are shown, 
constant star formation (solid line) and a declining r = 1 Gyr (dashed line). The amplitude of the difference between metallicities is related to the continuum 
color and the star formation history, where the absolute difference is largest for reddest colors and the shortest e-folding timescales. The dotted line shows zero 
change for reference. The histogram shows the observed colors of local blue sequence galaxies from the NFGS. 



The correlation of gas-phase reddening with the broadband U — V slope 

The effect of dust on the U-V colors of the stellar continuum (U - V)star depends on the adopted reddening law. Using the 
standard formulation (e.g., ICalzetti et al.ll200 ^. A\ = kx E(B - V) where A ^ is the dust attenuation of the stellar continuum at 
wavelength A, E{B - V) is the color excess, and k\ is the reddening curve, we have 

iU-V)ots= (U-V),,ar + iAu-Av)s,ar (Al) 

= (f/- VW + / (ku-kv) E{B-V)Hn (A2) 
where the color excess of the stellar continuum is linked to that of the nebular gas emission lines in H II regions via 

E{B-V)Hn = fE(B-V),„ar (A3) 

The ratio / represents a difference in the geometrical distribution of the dust. Estimates for local late type/irregular galaxies 
range from/ = 0.5 -0.1 (e.g., Kennicutt 1992, 1998; Calzetti et al. 1996; Bell & Kennicutt 2001). We leave it as a free parameter 
here and solve for it later. Note that we assume / does not depend on My, as is commonly done (e.g., Calzetti et al. 2000, 
Tremonti et al. 2004), but this is an important assumption which should be addressed in future studies. 

Differentiation of Eq. IA2l with respect to My yields the slope 



Sv{U-VU,s, = f (ku-kv) SvE(B-V)Hn 



(A4) 



where the slope 6vE(B - V)Hn follows directly from our Ha and Hj3 data and ku-kv = 
law. 
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The correlation of Ha equivalent width with the broadband U — V slope 

The dust-free line flux of Ha measures the instantaneous star formation rate (e.g. jKennicuttlfTgQStlKewley et alJl2002h . while 
the continuum flux measures the integrated star formation history. We can measure the dependence on broadband color of their 
ratio, the dust-free equivalent width, using the Ha-to-SFR conversion of Kennicutt (1998) 

SFR(MQyr-^) = 7.9 x 10-'*'LHa(ergi i"') (A5) 
and lBruzual & CharlotI (l2003h stellar population models to calculate the luminosity and colors. We then find 

logiWHa) = -0.83(±0.018)(t/-yW + 1.47(±0.08) (A6) 

which has deviations less than 5% over the range logiWna) =1-2 for all constant and declining star formation histories with 
r > 1 Gyr. We show in Fig.[T]that a linear approximation holds for a variety of star formation histories, and that the relation to 
color is tighter over the range log(W//Q) =1-2 than the relation to stellar population age or specific star formation rate. 

In the presence of dust, logW//Q will change if the stellar continuum and H II regions are obscured differently: 

logWffa = log(Wffa,ofo) + (AHaMII ' Aa6563)/2.5 (A7) 
= logiWHa.ohs) + 0.4 ^A6563 (1 - /) E{B-V)hII (A8) 

where WHa.obs is th e dust-reddened equivalent width and ^a6563 = 3.31 is the value of the reddening curve at A6563A for a 
ICalzetti et all (120001) reddening law. 

Combining Eqs A4, A6, and A7, and differentiating with respect to My gives us the slope 

-O.S36v{U-V)„ar = 6v \ogWHo.,obs + 0.4 kHo. (1 - /) 6vE(B-V)hii (A9) 

where we directly observe SvlogWHa.obs and 5vE{B-V)hii- 

The correlation of stellar metallicity with the broadband U — V slope 

We estimate the m etallicity Z through the reddening insensitive N2 = log(N 11 /Ha) indicator. Using the caUbration of 
iDenicolo et al] ( 120021) in terms of \og{0/H) 

12 + \og{OlH) = 9.12 (±0.05) + 0.73 (±0.10) A^2 (AlO) 

8.69 as the Solar abundance (lAUende Prieto et ani2001h . we find that the metallicity Z increases nearly 

Z/Zq =-0.22(±0.018)Mv- 3. 32(±0.25) (All) 
with a slope SyZ = -0.22 in Solar units. The metallicity ranges from ^ 0.4 Zq solar at My = -17 to ^ 1 .5 Z0 at My = -22. 

Unfortunately, the dependence of the broadband U -V slope on metallicity dv{U -V)z is a function of {U-V)siar rather 
than its derivative 6v(U -V)stai- This means we cannot relate it directly to the equations for dust and Whu and solve for it in 
a linear system in terms of SviU -V)siar- The solution is to take an iterative approach: we solve the system for the maximum 
(fixed) contribution of metallicity, recalculate Sv{U - V)z, and solve again with the updated value, iterating to convergence. 
This works well because the contribution of metallicity is never dominant (see e.g.. Fig. [T]). Even for a rather extreme 
r = 1 Gyr SFH and metallicities ranging from Z/Zq = 0.2-2.5 it never exceeds 6viU-V)z ~ -0.036 (or ^ 45% of the total 
observed U-V variation). The likely contribution is much smaller. In addition, if / > - meaning that there is at least some 
dust reddening of the stellar continuum - the effect of metallicity will decrease more as the absolute {U- V)siar will become bluer 



Putting it together 
We know that 

SviU-VUs = SviU-VUr + Sv(U-V)d,„ + SviU-V)z (A12) 

We form a linear system Ax = b from Eqs. 9 and from Eqs. 4,12 with SviU-V)srar and / as the unknown vector. The linear system 
is 

/-0.83 OAkHo.SvE(B-V)Hii \ ( SyiU -V)„ar\ ^ ( 5v\og{Wx{Ha)obs) + QAkHJvE{B-V)Hn 
\ 1 {kij-kv)5vE{B-V)Hn )\ f ) \ 5v{U -V)„b,-5v{U -V)z 

where we observe 



and with 12 + log(0///) = 
linearly with My 
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-0.078(±0.014) 


(A14) 


5vE{B-V)hu = 


-0.0060(±0.007) 


(A15) 




0.032(±0.016) 


(A16) 


SvZ = 


-0.22(±0.0018) 


(A17) 


We obtain the solution by solving the system once, using 


; the upper limit Sv{U - V)z ~ 


-0.036 for the contribution of 


metallicity, updating Sv{U - V)z, and iterating until converg 


ence. The solution converges 


quickly in four iterations with 



Sv{U-V)z = {-0.036,-0.013,-0.0068,-0.0046,-0.0044} for iteration {0, 1,2,3,4} respectively. The effect of metallicity van- 
ishes, because the best-fit {U-V)star is a constant « 0.25 over the whole My magnitude range (see Fig. Al). So we are confident 
to put in a more realistic initial Sv{U-V)z ~ -0.02 (SFH r = 5 Gyr and Z/Zq = 0.4- 1), and solve the system again to obtain the 
most likely answer 
The set of equations I A 1 3 l has the solution 



/ = 0.64(±0.14) (A18) 
6v{U-V)s,ar = -0.004(±0.015) (A19) 

Summarizing, the observed slope (5v(t/ -y)oij = -0.078 ±0.0014 has the following contributions from dust, stellar continuum, 
and metalhcity 

Sv(U-V)du,, = -0.073 ±0.016 (A20) 
Sv(U-V)s,ar = -0.004±0.015 (A21) 
Sv(U-V)z = -0.0014±0.003 (A22) 

(A23) 

We conclude that the local blue color-magnitude relation can be explained almost exclusively by a trend of increasing dust 
reddening with increasing absolute V-band luminosity. Surprisingly, the dust-free stellar continuum is nearly flat (U - V)siar = 
0.29-0.32 over the magnitude range My = 17-22. This result is not sensitive to choice of SFH and holds for most common dust 
laws. 
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